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SUPERSONIC FLOW INVESTIGATIONS WITH 
A “HYDRAULIC ANALOGY” WATER CHANNEL 
by 


JOSEPH BLACK, M.Sc., A.F.R.Ae.S., 
and 


O. P. MEDIRATTA, M.A., B.Sc. 
(Department of Aeronautical Engineering, University of Bristol) 


SUMMARY 


This paper describes the construction of a water channel in the University of 
Bristol for the investigation of the analogy between the two-dimensional flow of a 
gas and that of shallow water with a free surface. Both continuous and 
discontinuous flow were examined, with a view to determining the limitations of 
the analogy. 


Continuous “ shooting” water flow was found to be reasonably analogous with 
supersonic isentropic gas flow, a static depth of about half an inch appearing to be 


satisfactory with this particular channel. No independent check was made of the 
agreement, or otherwise, between “streaming” water flow and subsonic gas flow, 
since the method of checking used was the measurement of the angle of the waves 
formed on the water surface, and such waves exist only in “shooting” flow. 


Quantitative agreement between hydraulic jumps and shock waves was not 
obtained; this is in accordance with the theoretical predictions given in the text. 
It was found possible, however, to use hydraulic jumps to simulate qualitatively 
most of the features of gas flow with shock waves, such as flow in a Laval nozzle, 
or past a supersonic aerofoil; photographs of these jumps are included in the paper. 
The most valuable aspect of the water channel undoubtedly lies in the ease with 
which these phenomena may be demonstrated visually, and its construction for 
this purpose can be recommended. 


Further work of a quantitative nature, especially concerned with mixed flow 
regimes, is contemplated. 


1. INTRODUCTION 


The resemblance between the waves formed by a body moving at a supersonic 
speed through air, and the surface waves produced by a moving ship, or in water 
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flowing past bridge piers, was first commented on by Mach". It was not until 1920 
that the analogy* between two-dimensional gas flow and the flow of water with a 
free surface was first developed by Jouguet"’, and later by Riabouchinsky®’. They 
showed that the mathematical equations governing the two flows were of the same 
nature if the water were replaced by a fictitious “hydraulic gas” with a ratio of 
specific heats of 2.0. In 1937, Binnie and Hooker“? investigated the flow in a 
channel with a constriction. Preiswerk’’, in his application of the methods of 
gas dynamics to the accurate determination of the water flow in a converging- 
diverging channel, provided a very complete account of the analogy; since he was 
primarily interested in the water flow as such, and not in its possible application 
to gas flow, his work is of limited value for the purposes of this paper. Attempts 
to confirm Riabouchinsky’s work by quantitative measurements of lift and drag 
in streaming (subsonic) flow were made in 1939", and investigation of the flow 
through a cascade of turbine blades by this method was made in 1942’, 


The National Advisory Committee for Aeronautics have had a water channel 
in operation for some time, the first investigation reported’ being a comparison 
with the subsonic flow about a circular cylinder. A later paper’ reviews the 
practical details of the operation of this water channel, and gives some comments 
on the design of water channels in general. An examination of the corresponding 
flow of water and air in a number of cases has been made in Italy", and work has 
also been in progress in the California Institute of Technology; using a different type 
of apparatus, the emphasis there is on the relations between hydraulic jumps and 
shock waves"". 


In view of the comparative simplicity of the apparatus, and the promising results 
obtained by previous investigators it was decided to construct a small water channel, 
one of the objects being to determine the limitations of the method. 


Notation 
a_ velocity of sound 
velocity of surface wave propagation 
specific heat at constant pressure 
water depth 
Froude number, V /(gd)! 
Mach number, V/a 
pressure 


absolute temperature 


velocity component along x-axis 


*A brief account of the analogy is given in Appendix I. 
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v velocity component along y-axis 
V velocity 

y ratio of specific heats, C,/C, 

A wavelength 

p density 


surface tension 
Suffix 0 denotes reservoir conditions. 
Suffix 1 denotes conditions in front of the shock or jump. 
Suffix 2 denotes conditions behind the shock or jump. 


. DESCRIPTION OF THE APPARATUS 


4 
2.1. The water channel 

For ease of construction the channel was made inside a standard galvanised 
steel trough 5 feet by 18 inches, and 15 inches deep (Figs. 1 and 2). The bed of the 
channel consisted of a horizontal sheet of clear plate glass } inch thick, the water 
being led up to it over a convex dam of copper sheet. From the sump, the water 
was returned to the reservoir by a small centrifugal pump driven by a 4 h.p. D.C. 
motor. To lessen the surging and turbulence, the water entered the reservoir through 
a tube of fine mesh gauze which extended across the tank; this tube in turn was 
shielded by a curved sheet, and additional gauzes were inserted across the reservoir. 


Two methods were employed to accelerate the water over the bed; one was a 
converging-diverging nozzle, and the other a sluice gate. In the first, designed 
primarily for studying the flow through nozzles, there was a fixed contraction of 
varnished mahogany, which brought the width down from 18 inches to just over 
2 inches at the throat. Attached to the wooden contraction, and extending down- 
stream from the throat, were two vertical walls of 2-inch stainless steel strip; the 
nozzle profile over this section could be altered rapidly by push-pull rods pivoted 
to the strip. Movement of these rods was effected by means of screwed rods which 
passed through gimbals in the tank walls. 

This method had the disadvantage that the stream at exit from the nozzle had 
only a maximum width of about 4 inches, which was too narrow for the study of 
flow patterns around bodies in the stream. As an alternative, therefore, the fixed 
contraction together with the adjustable walls was removed, and a sluice gate, which 
could be raised or lowered to give a variable width slit at the bottom, was fitted 
across the upstream end of the bed (Fig. 1(5)). After passing under the sluice, the 
stream was contained between two vertical walls which ‘ran the length of the bed; 
these walls were diverged slightly to allow for boundary layer thickening. 


With both methods a simulation of back pressure could be made by raising an 
adjustable dam from the sump so that it projected above the downstream edge of 
the bed, and thus blocked the flow. 
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In order to diminish the effects of capillarity, which play no part in the analogy, 
the water in the tank had 0.05 per cent. Teepol added to it. This approximately 
halves the surface tension. 


2.2. Measuring apparatus 


Water depths at any point in the field of flow could be measured accurately to 
+0.002 inch with a micrometer depth gauge carrying a probe. This gauge was 
mounted on a block which could be traversed along a trolley by means of a 
micrometer-screwed rod (Fig. 2). The trolley in turn could travel the length of the 
tank; the position of the gauge in the flow field could be accurately determined 
from vernier scales on the trolley, and along the trolley guide rails. 


These guide rails of selected one-inch square brass were attached to, and ran 
horizontally the length of, the tank walls, one being mounted with a surface 
horizontal, and the other with surfaces at 45° (Fig. 1). A single point of support 
on one side of the trolley was provided by a small ball-bearing attached vertically 
to the trolley and running along the horizontal-surface guide rail. The other two 
points of support on the opposite side also provided the lateral restraint; they each 
consisted of a pair of small ball-bearings mounted opposite each other on the trolley 
at 45°, so that they ran along the top inclined surfaces of the other guide rail. A 
spacing of 6 inches between the two pairs was found to give sufficient restraint. 


Deflection of an electrostatic voltmeter, mounted in series with the survey 
probe, gave a positive indication of its contact with the water surface. This was of 
particular value upstream of the throat, and in regions of flow where the capillary 
waves formed by the probe tip on contact with the surface are difficult to detect 
visually. For the investigation of boundary layers and velocity distributions, pitot 
and static tubes made of hypodermic tubing, and mounted on a similar depth gauge, 
were used. Pressure measurements from these tubes were made on a Chattock 
three-liquid type micro-manometer*. 


2.3 Photographic technique 


Various lighting systems had to be used to suit differing types of flow. In 
general the flow was illuminated from directly below with a parallel beam shining 
through a large rectangular aperture cut in the bottom of the tank. For some cases 
the surface itself was photographed, while in others it was the shadow of the flow 
cast on to a ground glass screen above the surface which was recorded. 


For recording surface wave angles which were propagated along the length of 
the nozzle, a strip of 60 mm. recording paper was held in contact with the under 
surface of the glass bed, and a photo-flood bulb about 4 feet above the water 
surface was flashed on and off. 


*See Appendix II. 
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3. INVESTIGATION OF CONTINUOUS (SHOCKLESS) FLOW 


3.1 Analogy relations 


The analogy may be summarised as follows* : 


Two-Dimensional Gas Flow | Water Flow (Hydraulic Gas) 

Density ratio p/p, Water depth ratio d/d, 

Temperature ratio T/T, Water depth ratio d/d, 

Pressure ratio p/p, (Water depth ratio)’, (d/d, 

Velocity of sound, a,=(yp/p)! Wave Velocity, (gd)! 

Mach number M,=V/a Froude number F,= V / (gd)! = [2 (d, — d)/d}* 
Subsonic Flow Streaming Water 


Supersonic Flow Shooting Water 


From the expressions for analogous pressure and density relations, it will be 
seen that the ratio of specific heats in the “ hydraulic gas ” must be 2.0 (Equation (12), 
Appendix I), whereas in air, for example, it is 1.4. This difference does not limit 
the method. In the water channel experiments an attempt is made to get agreement 
between the measured values of the water depths and those calculated by gas- 
dynamic theory for the “hydraulic gas,” and not those calculated for a real gas. 
If such an agreement can be obtained, then by further calculation using the 
appropriate values of y the hydraulic gas results can be converted to those of 
a real gas. 


It is very important to understand clearly the significance of the Froude number 
as here defined. In general, the Froude number is the non-dimensional quantity 
V /(gly', where V is the velocity of flow of the liquid at a representative point, g 
is the acceleration due to gravity, and / is any representative length. The particular 
velocity and length which will be used in this work is the local velocity and the local 
depth, so that the Froude number throughout this paper refers to the 


quantity V /(gd}. 


First consider the velocity of flow V. From equation (6) in Appendix I, it is 
given by [2¢(d,—d)]', where d, is the depth of water in the reservoir and d is the 
local depth. This relation can be checked by direct measurement of the velocity 
with pitot-static tubes and comparison with the appropriate depth measurements. 


The quantity (gd)! is of a more complex nature. It is actually the velocity of 
propagation of gravity waves on the water surface for conditions where the depth 
d is small compared to the wavelength A (equation (14) of Appendix I); if these 
conditions do not exist the propagation velocity includes a further term, and 
(gd) has no longer this particular significance. 


*See Appendix I. 
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Fig. 4(a). 
Water flow along a converging-diverging nozzle (d, = 1.97 in.). 


In shallow water flow, therefore, where theoretically the depth d is small 
compared to the wavelength A, the Froude number as defined above is the ratio of 
the local velocity of flow to the local velocity of propagation of gravity waves on 
the water surface. As such it is thus directly comparable to the local Mach number 
in a gas flow, which is the ratio of the local velocity of flow to the local velocity of 
propagation of sound waves. It must be stressed that it is only the Froude number 
defined as V /(gd)' in the stipulated conditions, and no other, which can be regarded 
as comparable to the Mach number. 


In previous investigations: *’ the Froude number in the form [2 (d, —d)/d]! was 
used for the determination of the optimum depth of water to be employed. The 
depth d at the point where the Froude number F is unity will be % of the reservoir 
depth d,; if it is assumed that F is unity at the throat of the channel then the depth 
there should be 4 d,. Accordingly, the depths in the water channels were adjusted 
until they satisfied this assumed condition, and these depths were then quoted as the 
ones to be used. 

The objection to this procedure is that in itself it does not reveal whether the 
defined Froude number does in fact represent the ratio of flow velocity to wave 
propagation velocity. In other words, it does not check experimentally whether 
the flow conditions are such that the wave propagation velocity is accurately given 
by (gd)'. Such a check, however, can be made. 


A point source of disturbance in a shooting flow of water will propagate surface 
waves, and the angle which the envelope of these waves makes with the direction of 
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Fig. 4(b). 
Water flow along a converging-diverging nozzle (d,= 1.61 in.). 


flow will have the value sin-' (1/F) only if the depth conditions are as postulated 
in the analogy. Measurement of this angle thus provides a simple experimental 
check on the optimum depth of the water. 


3.2. Determination of the optimum depth 

The first step in this work was to measure the water depths at a number of 
stations along the nozzle, downstream of the throat, with the reservoir depth being 
held constant. A strip of Perspex, through which metal probes projected at right 
angles, was placed horizontally over the nozzle, so that the probes lay along its 
centre-line; the probes were then pushed down until their tips just touched the surface 
at approximately the same intervals as the original depth measurements. The 
surface waves were then recorded on recording paper held under the channel bed. 
Measurements of depths and photographs were repeated independently for each 
condition examined. In addition, pitot-static pressure traverses were made at the 
same locations in some cases. 

The wave angles on the recording paper (Fig. 3, see page 241) had to be carefully 
assessed. In order to attain the maximum reliability, the wave angles from the two 
sets of photographs were measured independently by two observers, before the 
Froude number had been calculated from the depth measurements. 

Typical results are shown in Figs. 4(a) and 4(b). It will be seen that at the 
downstream end of the nozzle the agreement between the Froude number as 
determined by the wave angle, and that from depth measurements, is reasonably 
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good considering the order of accuracy of the method. In fact, it appears to be as 
good as that obtained in a comparable experiment in air’, where the values of 
Mach numbers from wave measurements were thought to be good only to within 
+3 per cent. This agreement was maintained over the range of reservoir depths 
available with this channel, that is, 1.6 to 2.1 inches. Results obtained with different 
profiles also gave the same measure of agreem<at. 


The Froude number at the throat as determined from depth measurements, 
which naturally is the only means possible since no “ Mach” wave is formed there, 
was found to be unity, agreeing with the results of earlier investigators, that is, the 
depth in the throat was ¢ of that in the reservoir. 


Immediately downstream of the throat the values given by wave measurements 
are very scattered, but are all consistently higher than those from the depth 
measurements. This point requires further investigation, which is impossible with 
the existing channel; a channel is required which could produce a flow at outlet at 
least one inch deep, at the same values of depth ratios d/d, as in these cases. Results 
from this flow might then indicate whether the non-agreement arises from the 
increased depth. 


Part of the discrepancy may arise from deficiencies in technique. When the 
wave angles are large, as they are for small Froude numbers, for example, about 
40° at F=1.5, it becomes increasingly difficult to distinguish them from the curved 
capillary waves. In addition, the slope of the water surface, which may be seen 
in Fig. 4 to have become appreciable, would introduce a further error, since the 
recording method is essentially a shadowgraph of the surface. It is possible also 
that the expression for the wave propagation velocity, which is developed on the 
assumption that the wave is being propagated along a horizontal surface, requires 
modification when the wave is actually being propagated on a sloping surface, as 
in this case. 

Measurements of the velocity of flow with vertical pitot-static traverses indicated 
a uniform velocity profile above the boundary layer, the thickness of which at the 
end of the nozzle was of the order of 0.1 inch with a water depth of 0.5 inch. The 
velocity measured above the boundary layer agreed with that calculated from the 
depth measurements. 


These results would seem to indicate that for a water depth of about 4 inch 
the Froude number determined from the expression [2 (d,—d)/d]! may be taken 
as analogous to the Mach number in isentropic gas flow. This depth, however, 
should not be regarded as the only one to give good agreement; theoretical work 
still in progress appears to indicate that the flow need not necessarily be so shallow. 
The advantage of being able to use a greater depth is that the boundary layer is 
likely to be a proportionally smaller part of the measured depth. In designing a 
new water channel, therefore, a depth of about one inch should be allowed for, if 
possible; application of the technique already outlined will rapidly determine the 
optimum depth for the particular channel. 
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The N.A.C.A. investigation on flow past cylinders‘? was made at depths which 
were Claimed to be equivalent to subsonic Mach numbers; this implies an acceptance 
of the depth relation for the Froude number in a region where no independent check 
similar to that already described can be made, since wave forms of the type examined 
do not exist in streaming flow. One feature of their results is the different flow 
patterns created by the same body at the same “ Mach” number, for depths varying 
from 4 to 2 inches (Figs. 13 and 14 of Ref. 8). By comparing these patterns with 
that produced by a cylinder in a high speed wind tunnel the agreement is found 
to be best for a static water depth between 0.75 and 1.0 inch. Interpretation of flow 
patterns in streaming flow should be made with extreme care, however, since there 
is a type of wave formation which can occur in such conditions, and to which there 
is no analogy in air”*.’*’, These waves are being examined further. 


4. INVESTIGATION OF DISCONTINUOUS FLOW 


4.1. Shock waves and hydraulic jumps 


In the flow of either a gas or water, it is possible for a discontinuity to exist. 
With a gas, the phenomenon is the shock wave, passage through which causes an 
abrupt rise in the pressure, density, and temperature of the gas, and a decrease in the 
velocity. In the case of water flow, the discontinuity manifests itself as a sudden 
rise in the depth, which is termed the “ hydraulic jump”*’. A common example 
of such a jump is the circular one formed in the spread of a jet of water impacting 
on a horizontal surface (Fig. 5, see page 241)*. 


It is interesting to compare these phenomena, and hence examine the extent 
to which they may be regarded as analogous to each other. For simplicity, only 
jumps and waves normal to the flow will be examined. 


Water: The hydraulic jump can only occur in shooting water, the flow being 
reduced to streaming conditions behind it. Eddies are produced on the face of the 
jump (Fig. 6), and therefore there is a loss in what might be termed the “ mechanical ” 
energy of flow (the total energy of course remains constant, the eddying creating heat 
which cannot be re-utilised in liquid flow). 


The flow through the jump is governed by the following equations. 


Momentum : ‘ (2) 


*Another type of hydraulic jump, in which the rise occurs in an undulating, rather than an 
abrupt, manner, can also be formed. It is not discussed here, because it arises only if the 
depth ratio is below 2.0, and such conditions were not examined. In any case, no theoretical 
explanation of the undulating hydraulic jump has yet been given. 
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Fig. 6. 
Flow relations across a shock wave and a hydraulic jump. 


where / represents the eddying losses in the jump. 


The “ mechanical” energy of the flow is the sum of the pressure energy and 
kinetic energy, and is measured by the total head pressure represented by d,. 
Fig. 6 shows how the loss in total head Ad, through a jump increases with Froude 
number; for example, if the jump is formed at an F of 3.0 the loss in total head is 
about 1.3 times the depth of water in front of the jump. From these equations it is 
possible by expressing V in terms of d and the Froude number in front of the jump 
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to show that the depth ratio d,/d, is a function of the Froude number alone. The 
expression is 


d,/d,=4[((1+8 F?}- 1]. (4) 
The relation between d,/d, and F is plotted in Fig. 6. Measurements of d,/d, 


from the jumps formed in this channel are also plotted; the agreement is excellent. 
Similar results were obtained by Bakhmeteff"*’, and confirm equation (4). 


Gas: A shock wave can only be formed in a gas flowing at supersonic velocity. 
The equations governing its propagation are 


Continuity : p,V,=p,V;, 
Momentum: . . . . (6) 


Density relation (from Rankine-Hugoniot): "t-te (8) 


The gas flow has pressure and kinetic energy; in addition, it has a third form, 
which is its internal energy. Comparing the energy equations (equations (3) and 
(7)) for the jump and shock, it will be seen that in flow through a shock wave there 
are no losses of a like nature to the eddying loss involved. It is true that there is a 
loss in the total head pressure in going through the shock wave, and this might 
appear to be comparable to the same loss in the jump. This is not the case, however, 
because while the total head pressure d, is a measure of the energy of the water flow, 
it is the total head temperature 7, which measures the energy of the gas flow. The 
drop in total head pressure through the shock wave is only an indication of the 
re-distribution which is taking place between the kinetic, pressure and internal 
energies of the gas. 


Since the energy transport is so different in the two phenomena, it is to be 
expected that the relations derived and verified for continuous flow, which were 
based on analogous energy equations, will not apply to the discontinuous flow. It 
follows that if a hydraulic jump has formed in the water flow, depth measurements 
downstream of it cannot bear the same simple relation to the gas density as in the 
continuous case, and the analogy in that form breaks down. It is thus meaningless 
to put a value of 2.0 for y in the Rankine-Hugoniot relation, in order to obtain 
density ratios which will agree with those deduced from the measured depths 
through a jump. 


This lack of agreement may be seen in Fig. 6. If the analogy relations held, the 
line of d,/d, against F would also represent the density ratio p,/p, for a shock wave 
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formed in a “hydraulic gas,” where y was 2.0. The curves obtained by using a 
y of 2.0 in the shock wave relation are also plotted, and the divergence between it 
and the jump depth ratio is apparent. 


It should be noted that the values of depth ratio through the jump are for 
water to which no surface-tension reducing agents were added. Such agents are 
reported" to have a considerable influence on the form of the profile of the 
hydraulic jump. Ref. 11 also records a lack of agreement between hydraulic jumps 
and shock waves in air, in an investigation using an apparatus in which jumps were 
produced by a rapid discharge of water into a shallow layer of static water. 


4.2. Qualitative agreement between shocks and jump 

In spite of the quantitative difference between a shock wave and a hydraulic 
jump there is a remarkable agreement between their propagation and resulting 
wave patterns. Some examples of these are discussed. 


4.2.1. Flow through a converging-diverging nozzle 

The behaviour of a gas flow in a converging-diverging nozzle as the receiver to 
reservoir pressure ratio is steadily reduced, is well known*. When the receiver 
pressure ratio is only a little below that necessary to produce critical conditions at 
the throat, a normal shock wave forms across the nozzle slightly downstream of the 
throat; this wave travels down the nozzle as the receiver pressure drops, until it lies 
across the nozzle exit. 


With further reduction in receiver pressure the normal wave branches into two 
oblique shock waves which extend into the receiver. These are intersected by 
a small normal wave, but this disappears as the receiver pressure drops, until only 
the two oblique waves extending completely across the flow remain. When the 
design pressure ratio is reached, the flow at exit should be free from any shock 
waves, although it may still be intersected with Mach waves. 


With a decrease of receiver pressure below the design value, the gas at exit 
is flowing into a region at a lower pressure: consequently the jet expands in the 
familiar Prandtl-Meyer expansions, that is, through fan-shaped expansions radiating 
from the corners of the nozzle exit. 


In Fig. 7 (see page 242) may be seen the effect of decreasing back pressure on 
the formation of hydraulic jumps in a converging-diverging nozzle. The back 
pressure was simulated by raising the dam up from the sump. The progress of the 
jump is identical to that already described for the shock wave. In Fig. 7(b) the normal 
jump has moved downstream from Fig. 7(a), and in Figs. 7(c) and 7(d) the normal 
intersection between the oblique jumps is well defined. Flow out of the nozzle 
was not parallel, as it appears in the photographs; the straight edges are those of the 
illuminating slot below. 


*See, for example, Ref. 17. 
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Fig. 5. 
Circular hydraulic jump formed in the spread of a jet 
impacting on a horizontal surface. 


Fig. 11. 
Fig. 3. Channel blockage caused by putting the 
Surface wave angles as recorded model at such a large incidence that a 
on recording paper. second throat has been formed. 
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Fig. 7. 
Effect of decreasing back pressure on the flow of water in a converging-diverging nozzle. 
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Fig. 8. 
The hydraulic jump formed across “ shooting water by a 
second throat. 


Fig. 9. 
Hydraulic jumps formed on a double-wedge aerofoil at zero 
incidence in a shooting stream. 
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Fig. 10(a). 
The hydraulic jump formed ahead of a “* choked” duct in a “ shooting” stream. 


Fig. 10(b). 
The hydraulic jump system which results when the throat area of the duct in (a) 
is increased to allow the total mass flow to pass. 
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With the decrease to the back pressure of Fig. 7(e), the normal intersection has 
disappeared. For studying the flow patterns with further reduction of back pressure 
it was found that shadowgraphs were better than photographs of the water surface. 
To correlate the two, therefore, a shadowgraph of Fig. 7(e) taken with a circular- 
beam illumination is given in Fig. 7(f). 


The design back pressure is reached in Fig. 7(g), in which there is a parallel 
jet free from jumps. The Mach waves criss-crossing the flow are clearly visible. 


Finally, in Fig. 7(h), the expansion of the jet into a region of lower depth is 
seen. The spreading of the waves from the corners is very pronounced, and the 
general honeycomb pattern compares almost exactly with those sketched by Courant 
and Friedrichs"® for gas flow. There is a remarkable resemblance between the 
shadowgraph of the water “expansion ” and that obtained in an infinite fringe inter- 
ferogram of an expansion in air“*’. It is interesting to compare Figs. 7(f) and 7(h),' 
because although the extension of the waves appears identical in both, the former 
are clearly shock waves, since the waves do not intersect each other without mutual 
interference. With the expansion, however, the intersection is followed by a 
smooth and continuous deflection of the waves. 


4.2.2. Convergence of supersonic stream 


A similar effect to that of a high back pressure in the receiver, as already 
discussed, may be produced by making a supersonic flow converge by means of a 
second throat. This causes a shock wave to form. 


Again this phenomenon can be reproduced in water flow. Fig. 8 (see page 243) 
shows the normal jump formed across the flow ahead of the second throat. The wave 
formation from the probe tip touching the surface ahead of the wave shows that the 
water there is shooting, while the absence of such waves behind the jump confirms 
that the jump reduces the flow to streaming conditions. 


4.2.3. Double-wedge aerofoil 


The jump system around a double-wedge aerofoil at zero incidence is shown 
in Fig. 9 (page 243). Both the leading and trailing edge shock waves of gas flow are 
reproduced as jumps, and the Mach waves emanating from the nose surfaces and 
ending in the leading edge jumps can also be distinguished. 

When placed at an incidence greater than the semi-nose angle, the system of 
shock waves and expansions of gas flow was similarly reproduced, although it is 
not illustrated here. 


4.2.4. Supersonic diffusers 


An important practical problem is the diffusion of a supersonic gas stream 
entering a forward-facing intake. This diffusion, which accompanies the reduction 
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of the stream to subsonic conditions, depends upon the shock wave system existing 
at the inlet. Two distinct systems can be produced, depending on the amount of 
convergence of the duct. As a simple picture of what may happen, the duct may 
be regarded as a reversed Laval nozzle, and therefore there is an exact ratio between 
the inlet area and the area of the “throat” farther downstream, for any given inlet 
Mach number. 


If the duct converges to a throat area which is below the design value, it 
cannot swallow all the mass of air which the inlet is capable of passing. There is 
consequently a “ piling-up” of air in front of the inlet and some air will spill over. 
The shock wave formed in this case will be a bow wave ahead of the intake, thus 
resembling the wave formed ahead of a blunt-nosed body, although in the latter 
case the wave would be farther ahead. 


When the “throat” area of the duct is large enough, all the mass of air which 
the inlet can pass will flow through the duct. The result of this new flow condition 
is that the bow shock wave can move towards the inlet and become attached to 
the intake lips. The shock system in this case will consist of the oblique shock waves, 
both internal and external, from the lips. Depending on the Mach number and lip 
angle, the internal waves may have a normal wave intersection. 


Figures 10(a) and 10(b) (page 244) illustrate this phenomenon with hydraulic 
jumps. In the first case, the duct is over-converged and an intense jump forms a short 
distance ahead of the intake. With a larger throat area this jump has moved and 
attached itself to the lips; the internal oblique jumps from the lips coalesce 
with a small normal intersection and penetrate each other to reach the opposite 
walls of the duct. 


It should be noted that, while the norma! jump ahead of the duct in Fig. 10(a) 
must obviously reduce the flow to streaming conditions, the slight divergence of 
the duct from the throat to the trailing edge is sufficient to accelerate the flow to 
shooting conditions again. Thus the features of exit from a Laval nozzle are 
repeated, namely, two oblique jumps emanating from the trailing edges. Over the 
external surfaces in Fig. 10(b) the wave formation resembles that on the double- 
wedge aerofoil of Section 4.2.3. 


These wave patterns in Fig. 10(b) are of interest also as examples of the 
supersonic biplane as proposed by Busemann. 


4.2.5. Tunnel Blockage 


A problem which can arise in a two-dimensional test in a supersonic wind 
tunnel is illustrated in Fig. 11 (page 241). When a model is placed in the working 
section, it effectively divides the section into two channels. With a change of 
incidence, the “exit” area of the lower channel, formed between the model’s trailing 
edge and the nozzle wall, is reduced. This reduction, if sufficient to choke the 
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channel, will then cause an almost normal shock wave to form ahead of the “ inlet,” 
that is, from ahead of the leading edge across to the nozzle wall. It is thus comparable 
to the first case of the supersonic diffuser discussed in Section 4.2.4. 


In this photograph, the double-wedge aerofoil has been placed at a large 
incidence in the water channel, and immediately a normal jump has formed from 
slightly ahead of the leading edge across to the tunnel wall. 


The formation of small oblique roots to the normal jumps, which appear in 
this and other photographs, is probably due to boundary layer effects. 


The wave system in the water flow in this case is strikingly identical to that 
recently photographed in a supersonic wind tunnel*’; even the oblique root at the 
interaction of the shock or jump is similar in both cases. While such agreement may 
be largely fortuitous, it does suggest that further examination of the interaction of 
hydraulic jumps and boundary layers might be of value. 


5. CONCLUSIONS 

The water channel provides a convenient analogous flow to the isentropic flow 
of a gas. Since this type of flow is reasonably completely understood, 
there is little further knowledge to be gained about it from this apparatus. 


In the field where additional information would be welcome, that is, flow with 
shock waves, the analogy breaks down quantitatively, and it does not yet appear 
possible to establish any fundamental relation between the flow of gas through a 
shock wave, and that of water through a hydraulic jump. 


Fortunately the qualitative agreement between the phenomena is good, 
and the water channel is undoubtedly an excellent piece of demonstration apparatus, 
as the photographs included in the paper confirm. It might be of value, therefore, 
for indicating the shock wave pattern to be expected in any new system, for example, 
oblique shock entries for supersonic diffusers. 


There remain the serious restrictions that the method only gives an indication 
about conditions in two-dimensional flow, whereas most practical problems are 
three-dimensional, and that the boundary layer effects in the gas may differ widely 
from those in the water. 

Further work is contemplated with this water channel; it is hoped to examine 
the conditions of mixed streaming and shooting flow in order to determine whether 
transonic investigations are possible. 
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APPENDIX I 


The analogy between 
water flow with a free surface 
and two-dimensional compressible gas flow 


Consider water flowing along a horizontal plane bottom, bounded by vertical 
walls and with an upper free surface at constant pressure. 


A correspondence is established between the flows of water and gas by 
examining the continuity, energy, and wave propagation relations for each. 


J 


WATER GAS 
Fig. 12. 
(a) Continuity equation 


Water: Assume the vertical acceleration of the water flow to be negligible 
compared to that of gravity. Let u and v be the velocity components along the 
x and y co-ordinates, where x and y are measured in the horizontal plane along and 
normal to the direction of flow. The equation of continuity for two-dimensional 
flow is 


0 
3, (du)t jy ‘ ‘ (1) 


where d is the depth. 


Gas: For steady two-dimensional gas flow 


where p is the density. 


From (1) and (2) the depth in water flow corresponds to the density in the gas, 
that is, 
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(b) Energy equation 
Water: With frictionless flow, Bernoulli’s equation for the stream tube shown is 
gro, 
or »« @ 


The static pressure at any point in the flow depends only on its depth below the 
free surface, that is, 


P.=pg(d,—z,) and p=pg(d-z). (5) 
Thus from equations (4) and (5) 
The maximum velocity the flow can attain exists at the point where d is zero, 
therefore V max’ =2 gd, , 
hence (V =(d,—d)/d,. . (7) 
Gas: For gas flow the energy equation is 
Tal, Vn’. 
where T is the absolute temperature and C, is the specific heat at constant pressure. 
Thus . ‘ (8) 


From equations (6) and (8), or (7) and (9), it can be seen that the depth ratio 
in water corresponds to the temperature ratio in the gas, that is, 


T/T,=4a/d,. 


Equations (3) and (10) give two more relationships. Considering the gas only, 
it has an equation of state, p=pRT. Substituting the equivalent density and 
temperature relations, 


Thus the square of the depth ratio in the water flow is equivalent to the pressure 
ratio in the gas. 
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Ve 
Fig. 13. 


Since the isentropic pressure-density relation for a gas is p=K p’, 


P/Po=(P/ po), 


which from equations (3) and (11) can only be satisfied for a value of 


The water flow may thus be regarded, not as directly analogous to a real gas, but to a 
fictitious “hydraulic gas” (as Jouguet called it) which has a ratio of specific 
heats of 2.0. 


Velocity of wave propagation 

Water: In this case interest is centred in the velocity of propagation of gravity 
waves on the water surface, and not in pressure propagation through the water. 
In general, the expression for the velocity of propagation of surface waves is 


422) and 


where A is the wave length, p is the density, o is the surface tension, and d is the 
depth"*’. If the depth d is small compared to the wave length A, the velocity of 
propagation of the gravity waves becomes 


The capillary waves, represented by the second term, thus do not in this case play 
any part in the analogy. 
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Gas: The analogous wave in this case is a sound wave, which is propagated 
through the fluid with velocity a, where 


Froude number and Mach number 


The criterion for the dynamic similarity of the two wave propagations is the 
ratio of the free stream velocity to the wave propagation velocity. Thus in the case 
of water, from equations (6) and (14) 


V /c=V /(gd}' =[2 (d, —d)/d}'=F , the Froude number, . . (16) 
while in a gas, 
V/a=M, the Mach number. . 


For values of F or M less than one the wave can be propagated in all directions, 
its propagation velocity upstream of the source being V —a, and V +a downstream. 
With F or M greater than one, however, the wave front is the envelope of the circles, 
in the case of water, or spberes in a gas, the centres of which are being carried 
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Chattock micro-manometer for liquid pressures. 
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downstream with velocity V. The semi-nose angle of this triangular or conical 
envelope will be sin-' (1/F) or sin-! (1/M). 


Thus a water flow with Froude number less than one should correspond to 
subsonic gas flow, and F greater than one to supersonic gas flow, that is 


Mach number or 


Froude number Water Gas 
<i “ Streaming ” Subsonic 
>1 “ Shooting ” Supersonic 


APPENDIX II 


Chattock micro-manometer for liquid pressure measurements 


This instrument is a development of the more familiar water-oil Chattock 
micro-manometer for air pressures. The basic principle is that the displacement 
of a column of mercury in a glass U-tube, caused by the excess pressure in the 
left-hand bottle (Fig. 14), is compensated for by tilting the whole tube, so that the 
mercury column is again in equilibrium. An indication of the null point is given 
by the return of the magnified image of a water-oil meniscus to the intersection 
of cross-lines on a screen. To form this meniscus, oil is allowed to enter the right- 
hand bottle from the top until its separates the water in the inner tube above the 
mercury from that in the bottle. The meniscus observed is the convex one formed 
on the ground lip of the inner tube. 


Since the pressure difference must be equal to the height the right-hand bottle 
has been raised above the zero setting, it is necessary to measure accurately its 
vertical displacement. In order to do this, the U-tube is mounted in a frame 
pivoted at one end, and tilted by means of a vertical screwed rod, rotation of 
which is effected with a large disc. The height the bottle has been raised is thus 
easily determined from the amount of rotation of this disc, which has a scale 
engraved upon its circumference. Complete revolutions of the disc are indicated on 
the vertical scale at the right-hand end of the base. The number of turns can be 
converted to pressure difference by a conversion factor taking account of the distance 
between bottles, screw pitch, and relative density of water to mercury; it is thus 
possible to measure pressure differences to +0.005 inch water. 


REFERENCES 


1. Macu, E. (1887). Photography of Projectile Phenomena in Air. Sitzungberichte der Wiener 
Akademie, p. 164, Vol. 95, 1887. 


2. Joucuer, E. (1920). Some Problems in General Hydrodynamics. Journal de Mathé- 
matiques Pures et Appliquées (Series 8), Vol. 3.1., 1920. 


252 


= 


al 


10. 


HYDRAULIC ANALOGY FOR SUPERSONIC FLOW 


RIABOUCHINSKY, D. (1932). On the Hydraulic Analogy to Flow of a Compressible Fluid. 
Comptes Rendus de I’ Academie des Sciences, p. 998, Vol. 195, 1932. 

(1934). Some New Remarks on the Hydraulic Analogy. Comptes Rendus de I’ Academie 
des Sciences, p. 632, Vol, 199, 1934. , 


BINNIE, A. M. and Hooker, S. G. (1937). The Flow under Gravity of an Incompressible 
and Inviscid Fluid through a Constriction in a Horizontal Channel. Proc. Roy. Soc., 
p. 592, Vol. 159, 1937. 


PREISWERK, E. (1938). Application of the Methods of Gas Dynamics to Water Flows 
with a Free Surface. Mitteilungen der Institut fiir Aerodynamik, No. 7, E.T.H. Ziirich, 
1938. Translated as N.A.C.A. T.N. 934, 935, 1940. 


—— (1939). Some Applications of the Method of Hydraulic Analogy. Publications 
Scientifiques et Techniques du Ministére de I’ Air, pp. 43-77, No. 144, 1939. 


—— (1942). Brown-Boveri Revue, Vol. 29, No. 1-3, p. 77, Jan.-March 1942. 


Oruin, W. J., LinpNer, N. J. and Bitterty, J. G. (1947). Application of the Analogy 
between Water Flow with a Free Surface and Two-Dimensional Compressible Gas Flow. 
N.A.C.A. T.N. 1185, 1947, and also Report 875. 


MatTTHews, C. W. (1950). The Design, Operation and Uses of the Water Channel as an 
Instrument for the Investigation of Compressible Flow Phenomena. N.A.C.A. T.N. 2008, 
1950. 


TEOFILATO, S. (1949). A Contribution to the Study of Similarity between Hydro and Gas 
Dynamics. Monografie Scientifiche di Aeronautica (4) January, 1949. 

TEOFILATO, P. (1949). Monografie Scientifiche di Aeronautica (5) 1947, (9) 1948, (10) 1949. 
Gitmore, F. R., Presser, M. S. and Crosstey, H. E. (1950). The Analogy between 


Hydraulic Jumps in Liquids and Shock Waves in Gases. Journal of Applied Physics, 
pp. 243-249, Vol. 21, March 1950. 


BERSHADER, D. (1949). Interferometric Study of Supersonic Channel Flow. Review of 
Scientific Instruments, p. 266, Vol. 20, April 1949. 


HAvELock, T. H. (1908). Propagation of Groups of Waves in Dispersive Media. Proc. 
Roy. Soc, A. 81, pp. 398-430, 1908. 


Lams, H. (1932). Hydrodynamics. Cambridge University Press, 6th Ed., Chap. IX, 1932. 


BAKHMETEFF, B. A. and Matzke, A. (1935). The Hydraulic Jump in terms of Dynamic 
Similarity. Transactions of the American Society of Civil Engineers, pp. 630-680, Vol. 101, 
February 1935. 


Courant, R. and Friepricus, K. O. (1949). Supersonic Flow and Shock Waves. Inter- 
science Publishers, Ist Ed., pp. 387-392, 1949. 


Ferri, A. (1949). Elements of Aerodynamics of Supersonic Flows. Macmillan & Co., 
Chap. 8, 1949. 


LIEPMANN, H. W. and Bryson, A. E. (1950). Transonic Flow Past Wedge Sections. 
Journal of the Aeronautical Sciences, p. 750, Vol. 17, No. 12, December 1950. 


BARDSLEY, O. and Mair, W. A. (1951). The Interaction Between an Oblique Shock Wave 
and a Turbulent Boundary Layer. Phil. Mag., Series 7, January 1951, Figure 10, p. 33. 


253 


8. 
t 
9. 
11. 
|| 


CONTRACTING DUCTS OF FINITE LENGTH 
by 
L. G. WHITEHEAD, Ph.D., 


L. Y. WU, Ph.D. 
and 


M. H. L. WATERS, B.Sc. 
(Queen Mary College, London) 


SUMMARY 


A method of design is given for wind tunnel contractions for two-dimensional 
flow and for flow with axial symmetry. The two-dimensional designs are based 
on a boundary chosen in the hodograph plane for which the flow is found by the 
method of images. The three-dimensional method uses the velocity potential and 
the stream function of the two-dimensional flow as independent variables and the 
equation for the three-dimensional stream function is solved approximately. The 
accuracy of the approximate method is checked by comparison with a solution 
obtained by Southwell’s relaxation method. 


In both the two and the three-dimensional designs the curved wall is of finite 
length with parallel sections upstream and downstream. The effects of the parallel 
parts of the channel on the rise of pressure near the wall at the start of the contraction 
and on the velocity distribution across the working section can therefore be estimated. 


Typical designs are presented for contractions of moderate length which give 
small adverse pressure gradients on the walls and a uniform velocity distribution in 
the working section. 


INTRODUCTION 


The design of the contraction cone of a wind tunnel centres on the production 
of a steady stream with a uniform velocity distribution in the working section. The 
velocity distribution is associated directly with the shape of the contraction at its 
upstream end, while oscillations in the flow often arise from boundary layer 
separation near the start of the contraction caused by a too rapid change of section. 
In the present paper a method of design is given for contractions with circular 
cross-section and for two-dimensional flow at speeds for which compressibility effects 
are negligible. It enables good uniformity of velocity distribution at the exit to be 


achieved and boundary layer separation avoided in contractions of moderate length. . 


Frictionless flow is assumed so that for a real fluid the shapes should be modified by 
a slight displacement of the wall outwards to allow for the boundary layer thickness. 


Paper received March 1950. 
{The Aeronautical Quarterly, Vo!. II, February 1951] 
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The method has been applied initially to the design of a particular wind tunnet 
which has a contraction ratio of 7.25 and most of the examples in the paper relate 
to this design. Much larger contraction ratios, however, are well within its scope. 

It might appear at first sight that little risk of boundary layer separation would 
arise in a contraction since a large pressure drop occurs as a result of the reduction 
of the cross-sectional area of the stream. In short contractions, however, as found 
by Batchelor and Shaw", a rise of pressure frequently occurs near the walls at the 
start of the contraction, followed by a rapid fall of pressure to a value lower than that 
of the final flow. Thus there are two regions of rising pressure where separation of 
the flow from the walls can occur, at the start and at the end of the contraction. 

The two-dimensional flow problem is first considered and a series of contraction 
shapes giving a transition section of finite length from a parallel channel of given 
width to a narrower channel is developed from a simple boundary in the hodograph 
plane. Suitable choice of the parameters leads to contractions with small adverse 
pressure gradients and uniform velocity distribution in the working section. In 
the second part of the paper the flow with axial symmetry through these contraction 
shapes is calculated approximately. It is found that adverse pressure gradients are 
usually less than those for the two-dimensional flow and that favourable pressure 
gradients are increased. Boundary shapes which are satisfactory for two-dimensional 
channels therefore should have some extra margin of safety when used for tunnels 
of circular cross-section, in spite of the larger contraction ratio. The approximate 
solution depends on the close similarity which exists between the shapes of the 
streamlines for the two flows and becomes exact in certain simple cases. The 
accuracy of the approximation has been determined in a particular example by 
comparison with a numerically accurate solution obtained by Southwell’s 
relaxation method. 


Notation 
length AB in the hodograph plane { 
b length DE in the hodograph plane 
h resultant velocity of two-dimensional flow 
k two-dimensional contraction ratio (see Fig. 1) 
q.q; | complex variables in the hodograph plane 
R _ radius of the circular boundary in the hodograph plane 
u,v velocity components divided by U 
U __ velocity upstream of the contraction 
V __ resultant velocity of the three-dimensional flow 
W potential function for two-dimensional flow 
z=x+iy complex co-ordinates in the physical plane 
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a velocity potential for two-dimensional flow 
8 — stream function for two-dimensional flow 
y. 5 constants defined in terms of a, b and R 
#, ¥ constants defined in terms of y, 8 and W=2+i8 
@ velocity potential for three-dimensional flow 
wv Stokes’ stream function for three-dimensional flow 


SECTION I 


TWO-DIMENSIONAL CONTRACTIONS 
1. 

The approximate solutions of the three-dimensional flow equation described 
in Section II require a knowledge of the two-dimensional velocity potential and 
stream function for the same boundary shapes. Existing methods are not entirely 
satisfactory for this purpose. Lighthill'*) gives a convenient method for determining 
the boundary shape for a chosen velocity distribution along the wall, but the 
calculation of the shape of the intermediate streamlines and the velocity distribution 
across the channel is difficult. The method given by Cheers‘? leads to contractions 
of too great an overall length and in addition only gives a curve approaching the 
parallel section at either end asymptotically. A third method suggested by Hughes’, 
as giving the most rapid change of cross-section possible without reversal of pressure 
gradient, has a spiral discontinuity on the boundary so that it must be modified 
arbitrarily for practical applications. The present method avoids these defects. 


2. Flow in the hodograph plane 

The flow is first found in the hodograph plane where the co-ordinates are the 
velocity ratios u and —v. This choice is especially suitable in duct problems for 
the flow along a channel is represented in this plane by the flow from a source to a 
sink, and it is often possible to choose a simple boundary and determine the potential 
function by the method of images. Points at which the pressure gradient changes 
sign can also be found by inspection, since the radial distance from the origin 
represents the resultant velocity in the contraction cone. 


The co-ordinates in the physical plane, shown in Fig. 1, are x and y. The 
channel has an initial width of 2k, the velocity at this section is U and the two- 
dimensional contraction ratio is k. The channel wall in this z-plane and also in the 
hodograph plane is conveniently divided into four sections. The straight wall AB 
parallel to the axis, is followed by a curved section BC along which the vertical 
component of velocity increases to its maximum value at C. Along CD w decreases 
to zero and DE forms one side of the final parallel section. The corresponding 
configuration in the hodograph plane, defined by 
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z-plane. 


is shown in Fig. 2. W is the complex potential function in the z-plane and u and v 
are velocity ratios obtained by dividing the actual velocity components by U, the 
velocity at the beginning of the contraction. The upstream end of the contraction 
therefore corresponds to the point (1, 0) and the downstream end to the point 
(k.0). The area between the curve BCD and the real axis transforms into the upper 
half of the contraction, AE corresponding to the centre line. The configuration 
chosen gives adverse pressure gradients along AB and DE but these appear to be 
unavoidable in a contraction with initial and final parallel sections. If adverse 
pressure gradients are to be entirely avoided, the boundary in the hodograph plane 
must lie on, or between, the two dotted circles of radii 1 and k centred at O, and 
this can only occur if A coincides with B and D with E. The boundary in the 
physical plane then becomes a continuous curve and the parallel sections are 
eliminated. 

The shape of the boundary BCD must be selected to give a rapid contraction 
and this clearly requires as large a wall slope as possible near C. The simplest choice 
of a semi-circle for BCD is unsatisfactory as only a small vertical velocity component 
occurs at C. If, however, BCD is made a quadrant of a circle, as shown in Fig. 3, 
the maximum vertical velocity component is doubled and a contraction of sufficient 
rapidity results. The pressure along the curved wall then falls continuously from 
B to D. 

The flow in the g,-plane from the source at A to the sink at E can now be 
written down directly, six image sources and sinks being required to give the circular 
arc CD and the line BC as a streamline of the combined flow. The strengths of the 
sources and sinks are all equal and must be identical with the flow per unit time 
along the channel, i.e. 2Uk. It is convenient now to shift the origin in the hodograph 
plane to B and to change the scale so that the circle is of unit radius. A new variable 
q is chosen therefore and is defined by the relation 
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Fig. 2. 
Hodograph plane. 


where R, the radius of the circle in the g-plane, is given by 
R=k+a+b-1, 


and a and b are the lengths AB and DE, respectively. In terms of the new variable 
q the potential function then takes the form 


Uk (q° — a*/ R*) — R?/a*] 


The velocity potential and the stream function for the two-dimensional flow are 
denoted here by 2 and 8, respectively, to distinguish them from the corresponding 
functions ¢ and y for the axially symmetric flow. On the wall BC of the contraction, 
q is imaginary, so that £ is zero; on the centre line q is real and lies in the range 
a/R<q<(R-—b)/R, so that B= —Uk. 


3. Calculation of the boundary shape 

The flow in the z-plane is found by integrating the relation between dz/dq and 
q obtained by combining equations (1), (2) and (3). This can be done by re-arranging 
the equation and splitting into partial fractions but, since many terms result, the 
subsequent calculation of the shapes of individual streamlines is lengthy. 
Alternatively it is possible to make numerical integrations to determine the shape of 
any particular streamline directly from equation (3) and this procedure is usually 
more convenient. 
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The change in z along a streamline may be expressed directly in terms of the 
change in the velocity potential z as 


6z=62/(dW /dz), 


since 8 is constant. Substituting from (1) for dW/dz, separating into real and 
imaginary parts and integrating gives 


u and v are found from gq for equal intervals of z by re-arranging equation (3) in the 
form 
7W 


where y and 6 are defined by 


y= {R/(R-5)}?, 


The factor (eu _s) / be _ 1) is real on the outer wall and also on the centre 
line of the contraction, but complex for other streamlines. It is expressed generally 
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therefore as 2(u+iv), so that the solution of the quartic equation (5) becomes 


Equation (4) can now be integrated numerically. The only difficulty arises on the 
channel wall at C, where there is a rapid change in the slope of the wall. In this 
region smaller intervals of « are required than for the rest of the field and it is 
advantageous to start integrations from C as origin and proceed along the curved 
wall in both directions from this point. 


4. Simplified channels without parallel sections 


The elimination of the straight walls AB and DE simplifies the configuration in 
the g-plane considerably and the integration of the equations in complex form gives 
the shortest method of plotting the boundary. With a and b zero (3) becomes 


and hence dW = 2Uk [= 
dq = Lq q-l q+il 


Also combining (1) and (2) gives for this simplified case 
dW /dz=U [(k—-1) 


So that dz/dq takes the form 


This may be integrated directly giving 


z= =| klogq—log(1—q)- log(1 +4)+ {(k-1) 


(9) 


and it will be noted that the last two terms cancel in the special case k=2. The 
origin in the z-plane lies on the axis of symmetry because all the terms on the right 
hand side of equation (9) are real for q real and less than unity. The boundary of 
the contraction cone can now be plotted from (9) since q is known. For the 
intermediate streamlines equation (7) must be used in addition to determine values 
of g appropriate to the particular value of 8 required. 


A typical contraction of this simplified type is shown in Fig. 4 together with 
the resultant velocity A along the wall and the centre line. The two-dimensional 
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Long two-dimensional contraction without adverse 
pressure gradients. 


contraction ratio in this example is 2.69, corresponding to the value 7.25 in a 
three-dimensional flow. 


Figure 5 compares a duct having a contraction ratio of 4 with similar shapes 
given by Lighthill’’ and Cheers’. The present method is seen to give a more rapid 
change of cross-section near the centre of the range, but all three curves approach 
the limiting values slowly at the upstream and downstream ends. A further 
comparison of the three shapes is made in Fig. 6 which shows the maximum slope 
on the boundary for a range of contraction ratios. 


5. Examples with parallel walls upstream and downstream 


Examples giving a transition of finite length between parallel sections of 
differing width will now be considered. The length of the curved section BD, which 
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Comparison of contractions designed by different methods. 


1-6 
1-2 
oT 
0-8 
—— PRESENT METHOD 
0-4 ——— LIGHTHILL 
—-—- CHEERS 
0 
2 4 ‘ 6 8 10 
Fig. 6. 


Comparison of maximum wall slope for different methods. 


262 


4 
3 
y 

2 
0 

| 
|| 


CONTRACTING DUCTS OF FINITE LENGTH 


2:0 
— »b=0 
x x x b=0-16 
Uk 
1-0 
Fig. 7. 


Influence of a on the length BC. 


increases as a and b are reduced, is governed by the values chosen for the velocities 
at Band D. A rough measure of the length BC is given by the change in the velocity 
potential 2 between these points and this is given from (3) by the relationship 


(10) 


This change in 2 depends mainly on the value of a and is plotted, in Fig. 7, for a 
contraction ratio of 2.69. The small effect of varying b is shown by the difference 
between the full line for b=0 and the crosses corresponding to b=0.16. 


The corresponding relation for zp — 


_2Uk, 
°8 TR? + a") [R?—(R 


(11) 


is plotted in Fig. 8 and the main variable in this case is b. Values for higher 
contraction ratios can be found from these curves by increasing a and b in proportion 
to k-1. 


The choice of a suitable value for a must be settled as a compromise between 
the length of the curved wall and the adverse pressure gradient near B caused by 
the velocity falling to U(1—a). A typical contraction is shown in Fig. 9, in which 
the curved wall BD is plotted. The value of k is 2.69 and the. total change of the 
velocity potential between B and D has been taken at 4Uk. For this example a 
is equal to b and their value is 0.055, so that a rise of pressure of about 11 per cent. 
of the total head occurs on the straight wall near B. 
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Fig. 8. 
Influence of b on length C D. 


The value of b, in addition to varying the adverse pressure gradient beyond D, 
influences the rapidity with which a uniform flow is produced in the high speed 
section of the duct. In order to examine this effect, a further example has been 
considered which is similar to that shown in Fig. 9 except that b has been trebled. 
The change in 2 is reduced to 3.35Uk and the overall length of the contraction is 
reduced to 5.75, on the scale used in Fig. 9, due to a shortening of the curved wall 
CD. The difference between the centre line velocity at the start of the parallel section 
and the final value Uk amounts to 2.5 per cent. compared with one per cent. for the 
longer contraction. When the velocity distributions are compared on the basis of 
equal distances from the start of the contraction, however, the shorter contraction 
with the larger value of b is found to be superior. For values of x>6 the centre- 
line velocity differs from its final value by less than half the amount that it does for 
corresponding points on the curve shown in Fig. 9. It appears therefore that a 
gradual approach to the final parallel section is not advantageous, since it delays 
the rise in the centre-line velocity. A considerable increase of b is not permissible 
however, as the larger adverse pressure gradient on the wall would be likely to 
cause separation. 


SECTION II 
THREE-DIMENSIONAL FLOW 

The approximate method described in this section enables the three-dimensional 
flow to be calculated through contraction cones, of circular cross-section, identical 
in shape with the two-dimensional boundaries in Section I. The velocity potential 
a and the stream function £ of the two-dimensional flow are used as co-ordinates and 
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Fig. 9. 
Two-dimensional contraction of finite length. 


the method is basically similar to that employed by Fox and Southwell” for the 
calculation of the transonic flow through a convergent-divergent nozzle. 

With earlier methods such as those of Tsien’ and Szczeniowski'? it is not 
possible to investigate the effects of the parallel section upstream and downstream 
of the contractien. Thwaites’*’ method, in which contractions and expansions are 
repeated cyclically, appears to over-estimate the drop in velocity on the wall at 
the start of the contraction and to under-estimate the rapidity with which the velocity 
approaches uniformity at the high speed end. Contraction cones designed on this 
basis therefore will be unnecessarily long. 


2. Equations for the stream function 
The equation for Stokes’ stream function y~ for axially symmetric flow takes 


the form 
a (lay 


265 


L. G. WHITEHEAD, L. Y. WU AND M. H. L. WATERS 


where y is the height measured from the axis of symmetry. If the variables x and y 
are expressed in terms of = and 8, the equation readily transforms into 


a(1av) , 2 


aa) * ap) ~ °° (13) 


This form of the equation is more convenient than (12) when 2 and 8 are known 
for the boundary shape considered, because the boundary conditions are simplified 
to Y=constant for the two values of 8 appropriate to the wall and the centre line 
of the contraction. The wall and the centre line are streamlines for both flows, so 
that it seems reasonable to assume that the two sets of intermediate streamlines will 
not differ greatly in shape. The rate of change of ¥ with 2 therefore will be much 
smaller than the variation with 8. Following Fox and Southwell the a derivatives 
are neglected entirely as a first approximation, so reducing equation (13) to 


Applying a similar approximation in the calculation of the resultant velocity V 
for the flow, the exact expression 


h? av 
v= + Ga) 


reduces to V= P . FS) 
Comparing (14) and (15) it will be noted that 

0 

=0, 


hence to this order of approximation, V /h, the ratio of the resultant velocities for the 
two-dimensional and the three-dimensional flows, is a function of « only. Therefore 


— =f(z)=-= é (16) 
and integrating along an z=constant line leads to 

v=f (2) | ydB. 
If ¥ is chosen as zero along the centre line, its value on the outer wall will be 


Uk? /2 so that integrating across the channel yields the value of f(a) in the form 
0 


Uk? 
=1@| 


—Uk 
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Fig. 10. 
Contour for evaluation of f(z). 


Closer approximations to the flow may be found by substituting for the first term 
in equation (13) values appropriate to the first approximation and solving this 
modified form. This procedure has not been adopted because it has been found, 
by comparing a representative case with a numerically accurate solution obtained 
by Southwell’s relaxation method, that the first approximation has adequate accuracy. 


The direct numerical evaluation of the velocity ratio f(a) from (17) involves the 
calculation of values of y at a number of points along each ¢ line. It is often more 
convenient to express the integral { y df in terms of the values of § x da taken along 
the wall and the centre line. Consider first the integral 


(x +iy)(d2+idB), 


taken around the contour ABCD shown in Fig. 10, AB lying along the centre line of 
the contraction and CD along the wall. Since there are no singularities in the 
integrand within the area considered the integral is zero, so that its real part is 
zero. Thus 


(xdz+yd8)=0 , 


and hence, since £8 is constant along AB and CD and z is constant along BC and AD, 
B Cc 


yaB | xdz+ | xda. 
D 


A 


c 


ydB= 


B 


mony 


Only one integration across the stream therefore is required and this is conveniently 
made along the line <=0. The further integrations with respect to « make use ot 
data already found in determining the boundary shape and velocity distribution from 
the hodograph plane. 
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Fig. 11. 


Three-dimensional velocity distribution for long contraction. 


3. Examples of three-dimensional contractions 


Examples of three-dimensional velocity distributions calculated in this way are 
shown in Figs. 11 and 12. Fig. 11 shows the wall and centre-line velocities for the 
simple contraction without parallel sections illustrated in Fig. 4, and Fig. 12 gives 
the three-dimensional velocity distribution for the finite contraction shown in Fig. 9. 
In each case the contraction ratio now becomes k?=7.25. For the first example 
there is no adverse pressure gradient in the two-dimensional flow and none arises 
in the three-dimensional case; in fact moderate favourable pressure gradients replace 
the small favourable pressure gradients existing in the two-dimensional flow at each 
end of the figure. For the second example the adverse pressure gradients for the 
axi-symmetric flow are no greater than for the two-dimensional flow and the 
favourable pressure gradients are increased. 


4. Exact solutions 


It is of interest to note that the approximate solution of (13) becomes exact in 
special cases for which y may be expressed as a product of a function of « multiplied 
by a function of 8. A simple channe! of this type is the converging-diverging 
hyperbolic nozzle for which 


y=cosh zsinf. 


Other exact solutions of flow with axial symmetry depending on the existence of 
a simple relation between y, « and 8 are given by Milne-Thomson”?. 
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Three-dimensional velocity distribution for finite contraction. 


5. Numencal solution by relaxation method 


The relaxation method has been employed to provide a numerical check on the 
results shown in Fig. 11. The equations (12) and (13) are of the quasi-plane 
harmonic type considered by Southwell"® so that the method need not be described 
here in detail. Relaxation methods have previously been employed in the calculation 
of the flow in a wind tunnel contraction by Batchelor and Shaw who solved equation 
(12) for a particular boundary. In the present case equation (13) has been employed 
on account of the much simpler boundary conditions. The area over which the 
relaxation solution has to be determined in the z, 8 plane is a rectangle and all 
irregular stars are avoided. Equation (13) may be arranged in the form 


V? —3yV?y-yV? y=0, 


0? 
where V? denotes — + The residual F, at a point in a square mesh, is 


da? OB?" 
therefore given by 


4 


F= (yy) 3y, y+l2y,y,, 
1 


Me 


where the suffix 0 denotes the centre point of the star and the suffixes 1 to 4 the 
outer points. The corresponding relaxation pattern is shown in Fig. 13. Intervals 
of z and 8 of Uk/6 were chosen and the approximate solution was used to give the 
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(1, -3y,) 


(y,-3y,) ——— (lzy,-¥y) ——— (y,-3y) 


(1,-3y,) 


Fig. 13. 
Relaxation pattern for unit change of ¥. 


initial values of ¥. The small residuals were then quickly liquidated. The 
approximate and accurate values of ¥ are compared in Table I for the particular 
a line passing through C where the wall slope is a maximum and where the largest 
residuals occur. Very close agreement is shown, so that the estimated velocity 
distribution can be used with confidence. 


TABLE I 
Approximate value ay 0 29 114 256 453 703 1,000 
Accurate value ae O 28.7 1143 256.1 452.6 702.3 1,000 
CONCLUSIONS 


A method of two-dimensional design for wind tunnel contractions has been 
developed which gives a rapid change of cross-section with small adverse pressure 
gradients. The contractions are of finite length and the channel has parallel walls 
upstream and downstream of the curved section. An approximate method of 
satisfactory accuracy is given for finding the flow with axial symmetry through 
channels designed by the two-dimensional method. Boundary shapes satisfactory 
in two dimensions are also satisfactory for three-dimensional channels with axial 
symmetry. Estimation of the adverse pressure gradient at the start of the 
contractions and the uniformity of velocity distribution across the working section 
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require a method in which the parallel sections of the tunnel are represented. A too 
gradual transition to the parallel working section delays the development of a 
uniform velocity distribution across the working section and increases the length 
of the contraction unnecessarily. The use of a three-dimensional method which 
takes account of results found for two-dimensional flow leads to the design of 
satisfactory contractions of short overall length. 
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THE PROPERTIES OF CROSSED FLEXURE PIVOTS, 
AND THE INFLUENCE OF THE POINT AT WHICH 
THE STRIPS CROSS 


W. H. WITTRICK, M.A., A.F.R.Ae.S. 
(Senior Lecturer, Department of Aeronautics, University of Sydney) 


SUMMARY 


It is shown that the rotational stiffness of a crossed flexure pivot varies 
considerably when subjected to an applied force. The type of variation can be 
radically changed simply by moving the point at which the strips cross. The 
relation between torque and rotation for a given applied force is not exactly linear 
and the extent of the non-linearity is determined by taking into account the small 
movements of the centre of rotation of the pivot. Finally, for design purposes, an 
analysis of the maximum stresses in the strips is given. 


1. INTRODUCTION 


Crossed flexure pivots, of the type shown in Fig. 1, are extensively used in 
aeronautical research equipment to replace knife edges or ball bearings. Perhaps 
the most frequent use is for the fulcrum of the balance arm in a wind tunnel balance. 
Here, a null-point method of reading is almost invariably used, so that the prime 
considerations influencing the design of a flexure pivot for such a purpose are 
(a) that the applied loads should not cause instability of the pivot, and (b) that the 
rotational stiffness should be small enough to give adequate sensitivity. If changes 
in the applied load (e.g. the lift) should cause changes in the tension or compression 
in the strips of sufficient magnitude to give a noticeable change in rotational stiffness, 
the only effect would be a change in the sensitivity of the balance, with little or no 
change in the overall accuracy of the reading. 


There are many cases where changes in rotational stiffness affect directly the 
accuracy of the experimental readings. The present investigation, for example, 
arose primarily out of an experimental determination of the stability derivative /,. 
The technique employed consisted of forcing the model to perform rolling 
oscillations about its longitudinal axis, and the model had therefore to be mounted 
in such a way that these oscillations could take place. The obvious way to 
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implement this was by means of ball bearings. It was found, however, that the 
dust in the tunnel air made the ball bearings by no means frictionless, and in an 
attempt to overcome this they were replaced by flexure pivots, in which the strips 
were perpendicular to each other and crossed at their middle points. It was then 
found that the wide variations in the aerodynamic forces on the model gave large 
changes in the rotational stiffness of the pivots and their use was finally abandoned 
in favour of ball bearings. 

The behaviour of such symmetrical flexure pivots, in which the strips cross at 
their middle points, under a combination of vertical and horizontal loads, was 
discussed by the author in a previous paper“). It was assumed throughout that 
the rotation @ was small enough to be neglected in comparison with unity. With 
this approximation it was shown that the centre of rotation of the moving end of 
the pivot remained fixed at the point at which the undeflected strips crossed and, 
further, that the ratio between the shear force and tension in a strip was equal to 
6/2. The shear forces were neglected therefore and the tensions could then be 
obtained in terms of the applied load. These tensions were independent of the 
angle of rotation and the resulting relation between the rotation and the applied 
torque (about the centre of rotation) was linear. The analysis, confirmed by 
experiment, showed that considerable variations of rotational stiffness with applied 
force occur, and the pivot can become unstable, that is have a negative restoring 
torque, with the strips either in tension or compression, 


A recent paper by Haringx‘) discusses this identical problem, and derives 
results exactly similar to those mentioned above. He also considers the case of 
large rotations when the pivot is subjected to a pure torque only, with no applied 
forces. His results show that the ratio between the torque and rotation is not 
constant, but increases as the rotation increases. With the strips crossing at right 
angles the linear theory, assuming small deflections, underestimates the restoring 
torque by about 10 per cent. when the rotation is 45°. 


The present paper is concerned primarily with an investigation of the effect of 
changing the point at which the strips cross. It is shown that by adopting this 
simple expedient the variation of rotational stiffness with applied force can be 
altered radically. For example, by crossing the strips at one end, the pivot never 
becomes unstable with the strips in tension, although large variations of stiffness 
with applied force still occur. There is one pivot, however, in which the strips 
cross at 87.3 per cent. (=1/2+ ./5/6) of the distance along each strip, which appears 
to be superior in all respects. If the applied forces are not too great the centre of 
rotation remains fixed to a much higher order of accuracy than in other pivots. 
As a direct result of this the rotational stiffness remains sensibly constant over an 
appreciable range of vertical force. 

As a secondary consideration the linearity of the torque-rotation curve is 
investigated. It is shown that the departures from linearity are greater with, than 
without, a side force acting on the pivot. 
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The paper concludes with an analysis of the maximum stresses in the strips, 
for design purposes. 


length of each flexure strip 

total width of flexure strips on each side of the pivot 
thickness of each flexure strip 

inclination of each strip to the vertical 


a parameter whose value defines the point at which the undeflected 
strips cross (Fig. 1) 


bt*/12. Relevant moment of inertia of the cross section of the strips 
on each side of the pivot 


Young’s modulus of the material 

rotation of the moving end of the pivot 

lateral deflection of the moving end of a strip 
23/(@L)—1-p 

shortening of a strip due to change of slope 

e/(@? L). Given by equation (24) 

lateral deflection of a strip at a distance x from its fixed end 
x/L 

tension in a strip 

shear force in a strip 

[—PL?/(4ED}} 

value of y at which m, first occurs 

FL? /(EI) 

bending moment in a strip at a distance x from its fixed end 


moment of the terminal forces in a strip about the point at which 
the undeflected strips cross 


ML/(EI) 
stationary value of m. See Section 6 


vertical component of the forces applied to the moving end of the 
pivot, positive when inducing tension in the strips 


horizontal component of the forces applied to the moving end of 
the pivot 
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Fig. 1. 


Typical crossed flexure pivot. 


v VL? sec2/(El) 
HAL’ cosec 2/(El) 

T applied torque about the point at which the undeflected strips cross 
pB*/(B—tanh 

¢ (cothB—)+ ?6*/(8—tanh 

b 


coefficients in the relation between torque and rotation. See 
equation (27) 


K, value of K when V=H=0 


Suffixes A, B refer to the fixed and moving ends of a strip respectively. 


Suffixes 1, 2 refer to the strips on the two sides of the pivot. Strip No. | is that 
in which the horizontal force H is inducing tension. 


Suffix O refers to the point O at which the undeflected strips cross. 
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2. THE CENTRE OF ROTATION 


The pivot under consideration is of the type shown in Fig. 1. The strips are 
arranged symmetrically about both vertical centre lines, inclined at an angle < to 
the vertical and cross at the point O at a distance of 4(1+ )Z from the moving 
ends. Note that if p=O the strips cross at their middle points and if p= +1 they 
cross at the fixed or moving ends respectively. 


Consider now the pivot in a displaced position as shown in Fig. 2, the rotation 
of the moving end being @. The points B, and B, move to B,’ and B,’, the tangents 
B,’G and B,’G then making angles @ with the undeflected positions. The transverse 
deflections 5, and 5, are of the order of @L whilst the shortenings e, and e,, due to 
the change of slope of each strip, are of the order of 6?L. 


If the moving end is rigid, then B,B,=B,’B,’, giving 


(1+ p)? L? sin? z= { (1+ p)Lsin z—(6, —3,)cos z—(e, +e,)sin 2 } ? 


+ { (8, +5,)sin z—(e, —e,) cos 2 } *. 


If terms of the order of @ are neglected on the assumption that @ is small this 
expression reduces to 


2(1+p)L { 6, —5,)+(e,+e,)tanz } =(6,+8,)? cote. (I) 
Also, since the rotation of the moving end is @, 
(5, +5,) sin z—(e, —e,)cos z=(1+ p)Lsin sin 0. 
Neglecting terms of order 6° we have sin =6, and the last equation gives 
Squaring both sides of equation (1), it is seen that 
and similarly from equation (2) 
(4) 


Substituting (3) and (4) into the right hand side of (1) and neglecting terms of 
order 4° we obtain 


. (5) 
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Writing e=Aé?L, equations (2) and (5) when solved for 5, and 4, give 


6,=4(1+p)6L+ {4(1+p)tana+A, cot22—A, cosec 22 } 
§,=4(1+p)OL— {}(1+p)tan2+A, cot 22—A, cosec2z } 


It will be seen from equations (6) that if the terms in @ are neglected, then 


indicating that the point G coincides with O. Consequently, to the first order, the 
centre of rotation of the pivot is fixed. However, in order to deal with the non- 
linearity of the relation between torque and rotation it is necessary to retain the 
terms in 6’, thereby taking into account the movement of the centre of rotation. 


It should be noted that in the particular case of a pivot in which the strips cross 
at right angles (z=45°), equations (6) simplify to 


=4(1+p)OL+ (1+p)—-A, } PL 
(7) 


4 
(3 (1+p)—A, @L 


3. THE FORCES IN THE STRIPS 


Suppose that the moving end of the pivot is subjected to a system of forces 
whose resultant has vertical and horizontal components V and H through the point 
O, together with a torque T about this point, as shown in Fig. 3. The force V is 
positive when it is tending to induce tension in the strips but negative values must 
also be considered. The force H can be limited to positive values only, since the 
pivot is symmetrical about the vertical centre line and its behaviour as a whole is 
therefore independent of the direction of H. 


If P,, P, and F,, F, are the tensions and shear forces (referred to the undeflected 
directions) in the two strips, we have for equilibrium 


V=(P,+P,)cos 2+(F,—F,)sinz 


(8) 
H=(P, —P,)sinz—(F,+F,)cos 
It is convenient to express these equations non-dimensionally by writing 
VL? HL? 
v= sec, h= Ey COsee 2 
PL* FL? ©) 
48 EI f- EI 


where El (=Ebt*/12) is the total flexural rigidity of the strips on each side of 
the pivot. 
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Equations (8) then become 


+f.) cot . 


On solving for 8,? and 6, we obtain 


4B,?=4 (vu+h)+f, cot 22+ f, cosec 2a | 


(10) 
48,7 =4(v—h)-—f, cosec 2z—f, cot 2a 


In order to determine 8, and 8, completely, more information about the magnitudes 
of the shear forces is required. 


Let Mo, and Mo, be the moments about the point O of the forces in the strips 
at B,’ and B,’ respectively so that 
Mo, 
=M,,-4(1-p)F,L, 


with similar expressions for Mo,. Then the torque T is given simply by the equation 
T =Mo,+Mo:. . . (ly 


To proceed further it is now necessary to consider the bending of a single strip 
in detail, as shown in Fig. 4. The differential equation to be satisfied by the 
deflection y at a distance x from the fixed end A is 


Be 
Ela = Py+Msz Fx. 


On multiplying throughout by L’/(E/) and writing X¥=x/L, m,=M,L/(EID we 
obtain the equation 


48’y=L(m,—fX), 


where 8 and f are given by equations (9). 
The general solution of this equation is 


y=Q sinh 28X + R cosh 2BX + GX—-m,). . 2) 


The conditions of zero deflection and slope at the end A (X¥ =0). give 


_ -fL mL 
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Fig. 4. 
The bending of a single strip. 


Also, at the end B(X=1) we have y= and dy/dx=8, giving the two equations 


L 4? (cosh 28 — 1)— (sinh 28 — 2) 
f 

sinh 28 48 (cosh 28 — 1). 


On solving for m, and f we obtain 


8B? tanh 


BJ’ 
83 
R2 
_ 26 


It will be seen from equations (6) that € is of the order of 9. 


The bending moment M at any point in the strip is given non-dimensionally by 


_ML_1 dy 
x. 


By virtue of equations (12) and (13) this becomes 


m=m, cosh 2BX — sinh2BX. . . (18) 
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On using equations (14) and (15) this may be put in the form 


B nh BCI 2x) }. (19) 


m=6 B cosech cosh 8 (1 — 2X)+(p+ €) 


The bending moment M, at the end B, obtained by putting X¥ =1 in equation (19), 
is given non-dimensionally by 


B? tanh 8 } (20) 


ms=6 { Bcoth (p+ tanh B 


Also the moment M, of the terminal forces about the point O at which the 
undeflected strips cross, is given by 


Mo = MoL/(EI =m, —4(1—p)f. 
Using (14) and (15) this becomes 

where ¢=B (coth B—8)+p? 


and y is given by equation (16). 


Consider now the shortening e of the strip, due to the change of slope. To 
the accuracy required this is given by 


From equations (12) and (13) 


dy L 

9 26 

qX ag { 28m, sinh 28X — f (cosh 2BX — 1) } 
On substituting this expression into equation (23), performing the integration 

and substituting for m, and f from (14) and (15), we obtain after considerable 

simplification, 


| 


168 { B(3—coth? 8)+coth } +4(p+é)+ 


e=@L| 


B 


16(B—tanh By tanh 8 } |. 
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As pointed out previously, € is of order 6 and, for the accuracy required, can 
be neglected in the above expression. Hence, writing e=6*L, we obtain finally 


es a { 8(3—coth? 8)+coth } +4p + 


f 2 5 
* 76(6—tanh By | . (24) 


The forces in the strips have so far been derived in terms of 
B=[PL?/(4ED}. 
If the strip is in compression 8 becomes imaginary and for computation must be 
replaced by iy where 
y=[-—PL?/(4ED}}. 
All the functions of 8 contained in equations (14) to (24) remain real and can be 
obtained by the substitutions B=iy, sinh 8=isiny and cosh 8=cos y. 


It will be seen from equation (19) that the bending moment then contains 
a term proportional to cosecy which —>0°oO as y—>7, when buckling of the strip 
occurs. This imposes an upper limit on the value of y. 


The limiting values of the various bending moments and the shear force as 
8—>0 are required. These are as follows: — 


{ 1+3(p+)} 
{ 1—3(p+£)} 
f=60 (p+) » 
m=6 {14+3(p+é)(1—2X)} 
mMo=9 {1+3p(p+£)} 


Also the shortening A tends to the value 


PG. 
4. THE RELATION BETWEEN TORQUE AND ROTATION 


A study of the preceding equations shows that the determination of exact 
values for the tensions in the strips when the pivot is in a displaced position would 
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involve the solution of a set of simultaneous equations, each one being of a 
complicated transcendental nature. The torque T is clearly a non-linear function of 
the rotation @ for given values of V and H and the relationship will be represented 
here by a power series in 6, thus 


where the coefficients K, b, .. . are functions of v and A. 


Since terms in 6? have been retained, the equations so far derived are of a 
sufficiently high order of accuracy to determine the coefficients K and b in the 
series, but not the coefficients of the higher powers of 8. 


The magnitudes of K and b will be determined from the fact that, for given 
values of v and A, 


2 
K= EI\ G6), ‘ (28) 
1 L 


Consider first the strip tensions when @=0. From equation (15) it will be seen 
that the shear forces in the strips are then zero so that, by virtue of equations (10), 


B? (v+h) 
(30) 
B,?= 8 (v—h) 
and the two tensions are known. 
Further, from equations (11) and (28), 
_ fdmo, , dmo, 
Now, remembering that € is of the order of 6, equation (21) shows that 
dmo 
),..=9 
Hence from (31), 
K=9,+9., . . (32) 


where ¢, and 9, are the values of ¢ at 8=8, and 8, respectively, and 8,, 6, are 
given by equations (30). 
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BOTH STRIPS 
h 7 IN TENSION 


0 
ig. 5. 
The working region in the v-h plane. 
Let the value of K when V=H=0 be K,. Then since 8,=8,=0 from (30), 


Considering positive values of / only, it will be seen from equations (30) that 
if |v | >A both strips are in tension or both in compression according to the sign 
of v. If |v|<h/ one strip is in tension and one in compression. Further, one strip 
buckles in compression when v reaches the value (h—8z”). These considerations 
lead to Fig. 5 which shows the region in the v-h plane in which possible values of 
v and A must lie. It also shows whether the strips are in tension or compression. 
(Note that in Fig. 3, the strip No. 1 is the strip in which the side force H 
is inducing tension.) 


Equations (30) enable a chart to be constructed showing the values of 8,? and 
8,? for a range of values of v and h within the region of Fig. 5. The values of 
¢, and ¢, then follow from equation (22) and the stiffness K from (32). This has 
been done for three pivots in which p=0, +1 and + 5/3 (=0.745), and the results, 
given in the form K/K,, are shown graphically in Figs. 6, 7 and 8 respectively. 


In each of these diagrams the upper curve corresponds to H=0 so that the 
application of a side force always reduces the stiffness of a pivot. From Fig. 6 
it will be seen that a pivot in which the strips cross at their middle points (p=0) 
rapidly becomes unstable (i.e. the stiffness K becomes negative) when a vertical 
tensile force is applied. With a vertical compressive force the stiffness K rises to a 
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maximum of about 3.25K, (occurring at v = — 40) and then rapidly becomes negative 
as strip No. 2 approaches its buckling load. These characteristics are changed 
completely by crossing the strips at one end (p= +1), as will be seen from a 
consideration of Fig. 7. 


The application of a vertical tensile force then increases the stiffness whereas 
a compressive force decreases it. If h <62 approximately, the pivot never becomes 
unstable when v is positive. The intermediate value of p= + /5/3, the results for 
which are shown in Fig. 8, has been chosen because the value of 0K /dv at v=h=0 
is then zero. This pivot may be useful when a constant value for K is required 
over a small range of vertical force. 


Consider now the determination of the coefficient b in equation (27), giving 
a measure of the non-linearity of the torque-rotation curve. 


First differentiate equation (15) with respect to 9, remembering that € is of 
order 6, and then put 6=0, thus 


df 


(34) 


Now differentiating equations (10) with respect to 6, keeping v and / constant, and 
then putting 6=0, we obtain 


dé 


=- (scot cosec 22) 


B, 
d cot 22 cosec 22) 
(35) 


where y,, ¥, are the values of at and respectively, and 8, are given 
by equations (30). 


Next differentiate equation (21) twice with respect to 6, remembering again 
that € is of order 6, and then put 6=0, thus 


250-20} 


B tanh? 


(8 — tanh (37) 


where gus =coth — B (coth?8 + 1) + { 


Now, from equations (6) and (17), 


€,/06=4 (1+ )tan 2+2A, cot22z—2A, cosec 22 
(38) 


—£,/06=4(1+ ) tan «+2A, cot 22—2A, cosec 22 
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Hence, substituting (11) and (36) into (29) and using equations (35) and (38) we 
obtain finally 


b= { cos 22+¥3)— cos 22+y,)} (cosec 22)/4+ 
6) 2 


+(¥, —28,7) {4 (1+ p) tan «+ 2A, cot 2z—2A, cosec 22 } — 


—(¥, —28,7) { 4 (1+ tan 2+2A, cot2z—2A, cosec22}. (39) 


For given values of v and h, the values of 8, and £, are obtained from 
equations (30), then 9,’ and 9,’ from (37), ¥, and ¥, from (16), A, and A, from (24), 
and the value of b then follows from (39). 


If the side force H is zero, equations (30) give 8,=, so that ¢,’=9,’, ¥,=y, 
and A,=A, and equation (39) then gives b=0. This result could have been 
obtained simply by a consideration of the symmetry of the pivot since if the side 
force is zero the torque 7 must be an odd function of the rotation ¢. Consequently 
it is to be expected that the relation between the torque and rotation will be much 
more non-linear in the presence of a side force than with no such force acting. 


In the particular case of a pivot in which the strips cross at right angles (z= 45°) 
equation (39) reduces to 


(40) 


5. THE WORK DONE IN DISPLACING THE PIVOT 


In considering the variations of rotational stiffness with applied load it is 
instructive to investigate the work done by the vertical and horizontal loads V and 
H and by the torque T during a small rotation 6. For this purpose the horizontal 
and vertical movements of the point G (Fig. 2) are required in the directions 
of H and V respectively. Neglecting terms in 6° these may be shown to be 
4(A,—A,)@*L cosec a and | { }(1+p)—4(A, +A,) } sec respectively. 


It should be noticed that if h=0 so that 8,=8,=8 say and A,=A,=A 
say, the horizontal movement of G is zero and the vertical movement is 
4(1+p—4A)@Lsecz. If in addition v=0 so that 8=0, the shortening A is given 
by (26) and the vertical movement of G then becomes 


1 
2\ 92 
60 9p?) 6?L sec a. 


This is zero if p= + /5/3 so that in this case the point G coincides with O to a 
high order of accuracy. If |p| </5/3 the point G moves downwards (i.e. in the 
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direction of V) whereas if | p | >/5/3 the point G moves upwards. This explains 
the different types of stiffness variation shown in Figs. 6, 7 and 8. 


The work done by the forces V and H is given by 


Wy+Wu=@L[ V seca+4(A,—A,)H cosec 2]. 


Using equations (9) and (10) and neglecting terms in 6° this becomes 


6 
Wrt+Wu= 


L - +p- 4\,)+8,?(1+p—4,) } 


Further, neglecting terms in #°, the work done by the torque is given by 
W =(El6*/L) (40, +49). 
The strain energy U stored in a strip is 


1 
EI e 
U= 
0 
On substituting (19), performing the integration and neglecting terms in 6°, this 


gives after considerable simplification, 


42 


where ¢ and A are the functions of 8 given in equations (22) and (24). Hence 


giving a check on the preceding analysis. 


It should be emphasised that the point G is not the instantaneous centre of 
rotation of the displaced pivot but merely the displaced position of the initial 
centre. The actual instantaneous centre lies at a point whose horizontal 
and vertical distances from O in the directions of H and V are equal to 
{ QA,+A,)-—4(1+p)} seca and (A,—A,)@L cosec respectively, where terms of 
order 6? have been neglected. If p= + /5/3 and h=0 both these expressions are 
zero, so that for this particular pivot the distance between O and the instantaneous 
centre is of order 6’, and therefore the centre of rotation is fixed to a higher order 
of accuracy than in other pivots. 


6. THE GREATEST STRESSES IN THE STRIPS 


In designing a pivot the maximum stress in each strip, for given values of v, 
h and 6, is required. For this purpose the second order terms which have so far 
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been retained in the analysis can be discarded. This amounts to putting €=0 in 
equations (14) to (21), so that the terminal bending moments are then given by 


8? tanh } 

ma, n=O { B+p B—tanh B (41) 

and the bending moment m at any point in the strip becomes 

B*sechB . } 
6 Be +p 9 
m { cosech B cosh B (1 —2X)+5 tanh (42) 
The shear force f is given simply by 


The tensions in the strips may be assumed independent of the rotation for design 
purposes so that 6, and £, are then given by equations (30). 


In each strip the stress at any point is the algebraic sum of the uniform stress 
due to the tension and the stress due to the bending moment m. The maximum 
stress therefore occurs where the bending moment has its greatest numerical value. 
Because of the difference in behaviour of the hyperbolic functions which arise when 
8 is real and the trigonometric functions when f is imaginary, the tension and 
compression cases must be considered separately in determining the greatest 
bending moment. 


(a) The strip in tension (8 real) 
Consider first the case of p > 0. 


If the tension is zero (8=0) equations (25) with €=0 give 
m,=9(14+3p), mgp=6(1—3p), m=0(14+3p—6pX). 


Hence in this case | m,|>|m,|, and the variation of m with x is linear. The 
greatest bending moment is then equal to my. 


Further, since Bcoth8B>0 and f*tanhf/(G—tanhf)>0, it is clear from 
equation (41) that | m,| > |m,| for all real values of 8. This does not necessarily 
mean that m, is always the greatest bending moment, since for non-zero values 
of 8 the variation of m with x is not linear. There is therefore the possibility of a 
stationary value of m occurring inside the strip for some values of 8. If this does 
occur, however, it is clear from equation (18) that the magnitude of the stationary 
bending moment, m, say, is given by 
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This shows that | m, | > | m,| and it is concluded therefore that if p > 0 the greatest 
bending moment is m,, occurring at the fixed end, for all values of 8. An exactly 
similar argument shows that if p <0 the greatest bending moment is mz, occurring 
at the moving end. It will be seen from equation (41) that both possibilities may be 
combined in one expression for the greatest bending moment, max, by replacing 
+p by |p|, thus 


tanh 
{ 8 coth 8+ | p | 8—tanh B (44) 
(b) The strip in compression (y real) 
Substitution of 8=iy in equations (41) and (42) yields 
tan y 
ma, YOR (45) 
m=6 { ycosec y cos —2X)+p sin y (1 -2x)} (46) 


Consider first the case of p > 0. 


As before, if the compression is zero, | m,|>|msg{| and the variation of m 
with x is linear so that m, is the greatest bending moment. 


2 
Also, from (45), — =4p6 tany—y" 
Hence, if 0<y<7=/2, | m4 | >| me | 
and if /2<y<z, | |> | ma}. 


The possibility of a stationary value of m occurring at some point in the strip 
must now be investigated. If this happens its value m, may be obtained from either 
equation (18) or (46), thus 


2 J 2 2_y‘sec*y 
=6 {, cosec? y+ p (48) 
Using equations (45) and (48) we obtain, after some reduction, 
2 2 
m,? = mg, +0 —} . (49) 


Equations (47) and (49) show that if a stationary value does occur, it will be 
numerically greater than either of the terminal bending moments and therefore will 
be the greatest bending moment in the strip. 
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05 ie) 20 25 3-0 
lol 
Fig. 9. 
Variation of y,? with |p!. 


It remains to find at what value of y it first occurs. This is clearly the value 

at which dm/dX first becomes zero at one or other of the two ends of the strip. 

From (18), dm/dX=0 at X=0 when f=0. By virtue of equations (43) and 
(16) this gives 


If p > 0 this equation has a root, y, say, in the range 0< y<¢x/2 and for y>y}, 
a stationary value of m will occur inside the strip. 


Also, from equation (46), dm/dX =0 at X=1 when 


This equation has no root in the range O<y<rifO<p<l. If p>la 
root occurs in the range */2 << y <%, but this is clearly greater than y,, the root of 
equation (50), so that the stationary value of m, first occurs at y=7s. 


To summarise, it has been shown that if 0< y < y, the greatest bending moment 
is m, occurring at the fixed end, while if y> vy, the greatest bending moment is 
m, occurring at some point inside the strip. 

An exactly similar argument applies if »<0, in which case the greatest 
bending moment is mz if y is less than the lowest root of equation (51). If y is 
greater than this value the greatest bending moment is m,. 


As in the tension case the two possibilities of positive and negative p may be 
combined to give the following: 
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Fig. 10. 
Maximum bending moment in a strip. 


(i) if 0<y<js, where y, is the lowest root of the equation tan y=(1+| p|)y, 
the greatest bending moment occurs at one end of the strip and has 
the value 

|) tany 
4 YT tan ( 

(ii) if y> vy. the greatest bending moment occurs at some point inside the 

strip, its value being m, given by equation (48). 


The variation of y,” with |p| is shown graphically in Fig. 9. Also, the 
preceding results have been used to draw the curves of Fig. 10 which show the 
variation Of max/@ in a strip with the tension (87) for six positions of the point 
at which the strips cross, ranging from p=0 to |p|=1. In each curve the full line 
portion indicates that the maximum bending moment occurs at one end while over 
the dotted portion the maximum bending moment occurs somewhere inside the strip. 
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AN ANALYSIS OF THE LIFT ON STRAIGHT, YAWED 
AND SWEPT-BACK WINGS 
by 


J. LOCKWOOD TAYLOR, D.Sc., M.L.N.A. 
(formerly Professor of Aeronautics, The Technical University, Istanbul) 


SUMMARY 


As is well known, standard lifting-line theory cannot be applied to a wing 
with sweepback in order to determine the amount and distribution of the lift. 
Various methods, all fairly laborious, are available, and while they give results 
which agree fairly well with experiment, they do not show in detail the source of 
the various differences due to sweep. It is with the object of separating the several 
effects involved, which can be regarded as distinct for practical purposes, that the 
present method has been developed. Since one of the items, finite-chord loss, is 
shared with the straight wing, although not included in simple lifting-line theory, 
it is calculated first for the straight wing. The wing is represented by its median 
or camber surface and is assumed nearly flat. 


The various effects may be summarised as follows : — 


(i) Finite-chord effect; the reduction in lift for finite-chord (or aspect ratio) 
over and above that given by lifting-line theory. 


(ii) (a) Cosine effect at infinite aspect ratio; according to this the lift is reduced 
in the ratio cos A, where A is the angle of sweep. 
(b) Modification of (a) for finite aspect ratio applying approximately to 
either yawed or swept wings. 


(iii) The gain in lift at either wing tip, approximately equal to that for the rear 
tip of a yawed wing. 


(iv) The effect of the central kink in the plan form. 


Notation 
w  downwash velocity 
w, downwash velocity, lifting-line theory 
U _ aircraft velocity 
K circulation constant 
K’ _ strength of element of circulation 


Paper received January 1950. 
{The Aeronautical Quarterly, Vol. II, February 1951] 
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x co-ordinate measured downstream from origin of vortex 


y co-ordinate measured crosswind from origin of vortex 
X,, Yi» X2,¥_, cO-ordinates measured perpendicular and parallel to yaw or sweep 
respectively 
=x?+y?, + y’, + y’, =(y tanaA+ 2c} + y’ 
A angle of yaw or sweep 
semi-chord 
S$  semi-span 
A, aconstant 
6,% specified by scos@ or scos¢=distance from centre line to wing 
station 
integrals 
C, lift coefficient 
a angle of incidence 
a, lift-slope dC,,/dz2, infinite aspect ratio 
c.p. centre of pressure 
A.R. aspect ratio 


1. STRAIGHT WING FINITE-CHORD EFFECT 


The estimation of the effect of finite-chord (for a straight wing) is made by 
calculating the downwash at the } chord by an analytical method based on an 
exact solution for a uniformly-loaded rectangular wing. The derivation of this 
solution is being published elsewhere but will be summarised briefly. A single 
element of a wing vortex system of strength K’ comprising bound and trailing 
vortices, gives a downwash 


+= + 
U wy 


the origin of co-ordinates, x= y=0 being at the intersection of the two vortex lines; 
x is positive downstream and r?= x? + 


The second and third terms of the downwash represent the trailing vortex and 
the first the bound vortex. The combined vortex is anti-symmetrical in y and Tee- 
shaped, the vorticity at the x-axis being of opposite sign on the two sides of the y-axis. 
If to the first term is added a term 1/x, this annuls the vorticity at the axis of x 
for y negative and transforms the T vortex into a right angle ([). This can be 
regarded as the departure from two-dimensional flow round the aerofoil for 
negative values of y. Similarly, addition of —1/x gives the departure for positive y, 
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converting the T into *] (y positive to the right). For uniform chordwise loading 
over a wing of chord 2c 


= dx (y positive) 
K’ x=x+e 

=| 
_ K’ r,+y 
= (222) + y + y 
where r.= {(x+cP+y?}? 


The final term is the trailing-vortex downwash as used in lifting-line theory, 
the second term is the gradient of the trailing-vortex downwash, and the first term 
is due to the termination of the bound vortex. 


Numerical evaluation (see Appendix) shows that both the first two terms, 
representing the departure from lifting-line theory, can be represented with fair 
approximation, in terms suitable for integration, as follows, 


eqs 4x c 1 K’ 


2x c 1 K’ 
bound-vortex term ( y/e ) 
net downwash correction 2x £ ) K’ 

c y\1+2y/c/ 4 


The range of validity of the approximation to the net correction is x/y } 1.5 say, 
and for moderately small values of y, x }c. 


Only the net correction is required in the present application to a straight wing 
but the separate terms will be needed in dealing with yawed wings. Expressed as 
a fraction of the lifting-line downwash, the net correction is (x/c)/(1+2y/c) and 
if we assume this is independent of the chordwise distribution, provided x is 
measured from the centre of pressure (this has been checked by independent 
calculation for a parabolic chordwise loading (see Appendix) giving c.p.=0.25) 
and take x=c so as to give the downwash at the } chord (if the c.p. is as usual 
at 0.25), the integrals corresponding with 2c/[y(1+2y/c)] or 2c/y—4/(1+2y/e) 
must be evaluated. The first integral is the usual lifting-line integral, namely, if 
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K=4sU A, sin né@ in the usual notation, 


U cosé—cos¢ sin 
0 
The second integral is 
=- 6 10 
¢/2s + (cos 0—cos —¢]2s+ (cos 6—cos 
0 


the splitting into two parts being necessary in order to take care of the discontinuity 
at y=0 (or 6=¢@) already noted. The integrals can be evaluated without difficulty, 
and Table I gives, as an example, values for A.R.6, c/(2s)=1/ 12. 


The suffixes in Table I refer to n, the coefficient of @ in the series. Odd values 
only are required in symmetrical-load problems. 


Tas0—») =I't) (n odd) and Iqs-,) = —1%)(n even) gives the values for the other 


semi-span, required for the yaw case. 


The corrected downwash 


U 
TABLE I 

$ 0 30° 45° 60° 90° 
Re 0 1.037 1.055 1.030 1.015 
—2.4 — 2.37 — 1.223 — 1.108 —1.015 
I, 0 0.979 0.906 0.762 0.4155 
—1.6 — 1.415 — 0.216 0.020 0.4155 
0.812 — 0.157 — 0.448 
tT. — — 0.631 — 0.58 0.448 
I, 0 0.569 0.069 — 0.356 — 0.34 
—0.711 — 0.062 0.581 0.381 — 0.34 
I, — 0.29 — — 0.454 0.294 
— 0.304 — — 0.135 — 0.294 
I, 0 —0.03 — 0.435 —0.159 0.29 


— 0.316 0.50 — 0.05 —0.39 0.29 
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Hence equations corresponding with those of Glauert for the lifting-line can be 
written down by virtue of the relation 
w 


For the particular case of the rectangular wing of A.R.6, taking a,=27, we find, 
taking three terms of the series and using the 0, 45°, 90° stations, 


A,=0.226 A,=0.014 A,;= —0.0018 (a=1) 
or taking 30°, 60°, 90°, 
A,=0.2155 A,=0.006 A,;= —0.0023. 


It seems reasonable to take the mean values 


A,=0.221 A,=0.010 A,= — 0.002. 


The convergence is very satisfactory and considerably better than for lifting- 
line theory. The lift-slope is *sA,/c=0.663 x 2x, as compared with Glauert’s value 
0.72a,, a decrease of nearly 8 per cent. Figures in Ref. 1 from various sources 
range from 5 to 8 per cent. for similar aspect ratios, and there is experimental 
confirmation from two-dimensional tunnel results. 


It is of interest to note that the loading is nearly elliptic, even for the rectangular 
wing when the corrected downwash at the } chord is taken. The c.p. shift can be 
obtained by comparing the downwash w, at } chord (or assumed c.p.) and at 
} chord, and hence evaluating the induced camber. It is found that the shift is 
small except at the tip, the mean being 0.01 chord giving c.p.=0.24 corrected. 
The induced drag, since it depends on the vortex wake at a distance (i.e. K) will 
be very nearly as for elliptic loading in spite of the appreciable increase in incidence 
for a given lift. 


2. YAW EFFECT 

(a) The reduction of lift in the ratio cos A for an infinite-span constant-section 
yawed wing follows from resolution of the velocity U into components U cosA and 
U sind normal and parallel to the wing generators. The incidence, camber ratio 
and zero-lift angle are all increased in the ratio secA for a normal section as 
compared with one parallel to U so that the lift is reduced as cosA at any given 
fore-and-aft incidence. A similar result is obtained by considering the downwash, 
remote from the trailing vortex, due to a bound vortex strength K cosA distributed 
over the normal chord 2ccosA. The vertical velocity component at any point in 
the chord is as for vorticity K distributed over chord 2c, as required. 


(b) At finite aspect ratio there will be a correction to the cosA factor. This 
has been estimated for a 30° yawed uniform-chord untwisted wing of aspect ratio 6 
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(the wing tips are streamwise) by substituting (e—w/U)cosA for 2—w/U in the 
fundamental equation and allowing for sweep in the calculation of w for given K 
as follows. Using the same approximation as before for the trailing vortex field, 
but writing x/c=1-—(y tan A)/c to allow for sweep or yaw angle A, 


w_ K’ _ ytana 2 
A 
y 
x, 
x 
Fig. 1. 
Yawed axes. 


gives the downwash due to a trailing vortex element K’ at c (half-chord) downstream 
from the centre of pressure line, 


(¢ - - 
or (y positive) 
_ ytand (¢ 2 
Similarly ) 73) (y negative) 

_K 
\y -14+2y/ce7° 


The bound-vortex terms require separate treatment. Using co-ordinates y,, x, 
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parallel and perpendicular respectively to the sweep and therefore connected with 
y. x by the relations 


y,=(—xsinA+ycosA), x,=(xcosA+ysinA) (Fig. 1), 
the uniform-chordwise-load solution for the yawed wing is 
x+e 
K’ (2 — [ sec A log (r+ y,)+ tan A log (r-+2) | 


=) 


x=x+e 


x=x-c 


=— 
4x sec A log + tan A log "+x+e 


At the } chord line, x=c-—ytanA is to be substituted; the resulting expression is 
given in the Appendix. 


Numerical evaluation (see Appendix) gives as a working approximation 


4m y 


= /( )seca=- 2 ) seca (y positive) 


1+2y/ce 1+2y/c 


i.e. as for the unyawed wing but with the additional factor secA; there is a 
corresponding equation for y negative. 


Collecting results for both bound and trailing vortex elements 


w_K’ (2—secA)c 2(2+tanA—secd) rae 
{ 1+2y/e } (y positive) 
we K’, (2- _ 2(2+tandA—secd) 
2{ (y negative). 


The integrals to be evaluated are as before, apart from the constant factors depending 
on A, the angle of yaw, namely 


= =(3— sec A) — (2+ tan A - sec (2—tanA—secaA) SI’. 


Even values of n now occur in addition to odd values. Taking coefficients A, to A; 
inclusive and satisfying the downwash equation at ¢=30°, 60°, 90°, 120°, 
150° we get (for 30° yaw) A,=0.1969, A,=0.0199, A,=0.0119,. A,=0.0061, 
A,= —0.0005. 
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Fig. 2. 
Central kink axes. 


Comparing these values with those for the unyawed case, for which A,=0.221, 
A,=0.010, A,= —0.002, the main difference in the symmetrical terms is .a reduction 
in A, of about 11 per cent. to 0.89 of its original value, compared with 0.866 for 
cosA. The difference from infinite aspect ratio therefore is fairly small, and it 
seems reasonable to apply the result also to the swept wing of similar aspect ratio. 
The factor 0.89 compares with 0.85 given by Falkner“) for an elliptic wing of 
aspect ratio 6.36 at 30° yaw. 


3. TIP EFFECT 


The results just calculated for a yawed straight wing include the local change 
of lift at the tips, namely decrease for the leading tip and increase for the rear tip. 
The fractional increase of lift on the rear semi-span is 


[ (2/3) A, — (4/15) A,]=0.076. 


That this is essentially a local effect can be shown by plotting the lift-grading 
(Fig. 3) and by evaluating the downwash due to the A,, A, terms from the other 
semi-span, which is found to be negligible. Since this will apply equally to the 
swept as well as to the yawed wing it may be concluded that a similar increase of lift 
of 7.5 per cent. localised at the tip, will apply to each tip of a swept wing. The 
A, term also gives a value for the /, derivative (rolling moment due to sideslip), 
namely [ (=s/4c) A,/(=/6)]A=0.18A for the yawed wing. This is in fair agreement 
with other data (Ref. 3, p. 227) making allowance for some departure from linearity 
at the extreme sideslip angle of 30°. 


4. KINK EFFECT 


The central kink effect, again based on the exact result for uniform chordwise 
loading, may be evaluated by somewhat similar means. Considering a kinked 
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K 0-02 


Fig. 3. 
Local increase of lift at rear tip of 30° yawed wing, A.R.6. 


straight-line vortex and taking axes (x,, y,) (x;, y.) with the y’s along the respective 
vortex lines and the origin at their intersection so that 

y,=xsinAd+ycosaA 

x, =xcosA—ysinA 

y,=xsinA—ycosA 

xX,=xcosA+ysina, 
(x, y are measured downwind and crosswind respectively (Fig. 2)) 


then 
@eN\rx, 


gives the downwash measured from the two-dimensional flow for the port half-wing, 
represented by 2/x, as datum. Integrating as before 


x+e 
K 

|ax- ( ) + tan A log (r+2)| 

x=x-c 
K { sec log 27 + (c/y)tand) tan A+ 2c 

4 r+ytanaA+ ic 

(r’ sin A+ y tan A + 2c) (sec A + tan A) 

where r= {(ytand+2cP+y? } 
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Fig. 4. 
Loss of lift due to central kink (30° sweep). 


Numerical evaluation gives the following results 


AR. 2 4 6 
0.40 0.34 0.28 
0.10 0.085 0.07 


2zw/U is equal to the reduction (8C,) of C,, due to the kink. Since the actual 
C, at the centre line of the wing is 4K, the fractional loss of lift, averaged over the 
span, is as given. The losses referred to the mean wing C, are estimated as being 
20 per cent. greater. Fig. 4 shows the spanwise distribution of 8C,/C, 
due to the kink. 


Summing the various lift changes for a swept wing of A.R.6 compared with 
a straight wing 


Due to cosine effect (Section 2(a)) —11 per cent. 
Due to wing tip effect (Section 3) + 7.5 per cent. 
Due to central kink (Section 4) — 8.5 per cent. 


—12_ per cent. 


8C, /C, 

net 
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Ref. 2 gives a lift loss of 14 per cent. for an elliptic wing of similar A.R. and there 
is some experimental confirmation from Ref. 4, pp. 394-5. If comparison is being 
made with straight-wing lifting-line theory the finite-chord loss (i) of about 8 per 
cent, must be added, giving 20 per cent. in all. 
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Comparison of exact and approximate values of downwash 


correction 


(1) Straight wing 


Trailing-vortex correction 


Exact 


Bound-vortex correction log ( 


Cc 
Approximation to net correction, 2 <( 


x/c (measured from c.p.}= 0.5 


y/c=0.5 
y/e=1 | 
y/e=2 


yle=3 


trailing vortex 
bound vortex 
net correction 
approximation 


trailing vortex 
bound vortex 
net correction 
approximation 


trailing vortex 
bound vortex 
net correction 
approximation 


trailing vortex 
bound vortex 
net correction 
approximation 


1.74 
0.56 

1.18 (1.4) 
1.0 


0.68 
0.28 
0.40 (0.43) 
0.33 


0.22 
0.10 
0.12 (0.12) 
0.10 


0.103 
0.05 

0.053 
0.047 
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uniform chordwise loading. 


1.0 


3.12 

0.94 

2.18 (2.32) 
2.0 


1.24 
0.47 
0.77 (0.79) 
0.67 


0.415 

0.19 

0.225 (0.23) 
0.20 


0.20 

0.093 
0.107 
0.094 


1.5 


1.57 
0.56 
1.01 
1.0 


Figures in brackets are for parabolic chordwise loading. 


2.0 


0.68 
0.28 
0.40 
0.40 


3.0 


0.47 
0.19 
0.28 
0.29 


| 
| 
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(%) Yawed or swept wing 
The trailing-vortex correction has been taken as twice the approximation given 
for the net correction, which, as will be seen, limits the angle of sweep to 


rather less than 45°, x/y=1 approximately. The bound-vortex correction has been 


taken as secA *( 5), the values of which are compared with the exact 
y \l+2y/c 


expression, 


2c—y A 2) — inA A 
sec A log { V[(2c—y tan 2esin +y sec 


+tan tog { y tanA) + y?]+ 2c—y tan 


y (sec A — tan A) 
y/c= 0.5 1.0 2.0 4.0 

Exact 1.41 0.65 0.22 0.045 

Approximate 2.0 0.67 0.20 0.057 


REFERENCES 


1. FALKNER, V. M. (1943). Aerodynamic Loading on Surfaces of Any Shape. R. & M. 1910, 
August 1943. 


Couen, D. (1945). The Twist of a Surface for a Given Lift. N.A.C.A. Report 826, 1945. 


CAMPBELL, J. P. (1948). Factors Affecting Lateral Stability. N.A.C.A. University Conference 
on Aerodynamics, June 1948. 


RICHARDS, E. J. (1947). Practical Design Problems Arising from Sweepback. Aeronautical 
Conference London, 1947. Royal Aeronautical Society, 1948. 


304 


é 
! 
‘ 
q 


NOTICE TO CONTRIBUTORS 


All communications should be addressed to the Secretary, Royal Aeronautical 
Society, 4 Hamilton Place, London, W.1. 


Papers for consideration should be typewritten with double spacing and be 
accompanied by a summary of not more than 250-300 words. 


Only simple symbols and formulae should be typewritten. All others should 
be written clearly in ink. Greek letters should be indicated and a clear distinction 
drawn between capital and small letters. Particular attention should be paid 
to the use of alpha and a, kappa and k, mu and wu, nu and v, eta and n. The 
letter / must be distinguished from the numeral 1, and the letter O from the 
zero 0. 

The use of dots, bars, and so on, over letters or the use of dots as multi- 
plication signs and bars for brackets should be avoided as far as possible. 

Suffixes and prefixes must be clearly indicated and complicated suffixes 
and prefixes should be avoided. The square root sign / should be replaced 
by the index 4. 


A list of all symbols used must be given at the end of the paper whether they are 
standard or not. 


Clear diagrams or sketches must be submitted for redrawing by the Society’s 
draughtsman. 

Photographs must not be less than half plate size and must be clear black 
and white glossy prints. 


References should be given in the form:— 
Warp, G. N. (1949). Calculation of Downwash Behind a Supersonic Wing. 
The Aeronautical Quarterly, Vol. 1, May 1949, pp. 35-38. 


The copyright of every paper printed in The Aeronautical Quarterly is the 
property of The Royal Aeronautical Society. Permission to reprint or to use 
any paper will not be unreasonably refused. 


It is a condition of publication that the author shall have obtained consent where 
necessary to use any material in his paper which is copyright or the property 
of any other person or his employers. 


Authors will be entitled to 25 reprints of their papers free of charge. 


The Aeronautical Quarterly is published at the Offices of The Royal Aeronautical Society, 
4 Hamilton Place, Piccadilly, London, W.1, at 10s. 0d. net for a single number to non-members 
and 7s. 6d. net for a single number to members of The Royal Aeronautical Society. Instructions 
for Contributors are given above. None of the papers must be taken as expressing the 
opinion of the Council of The Royal Aeronautical Society unless such is definitely stated to 
be the case. All communications for publication or advertisements in The Quarterly should 
be addressed to—The Secretary, J. Laurence Pritchard, 4 Hamilton Place, W.1. 


Telephone: Grosvenor 3515 (5 lines). Telegraphic Address: Didaskales, Audley, London. 


THE AERONAUTICAL QUARTERLY 


VOLUME II 
PART I, MAY 1950 
Two-Dimensional Theory of Stiffened Plates . 
The Whirling of a Spinning Top 


The Possibility of the Determination of Rate of Climb 
from Acceleration Measurements in Level Flight ... 


A Note on Subsonic Aerofoil Theory 
Theory of an Oscillating Supersonic Aerofoil 


The Linearised Theory of Conical Fields in Supersonic 
Flow, with Applications to Plane Aerofoils es 


PART II, AUGUST 1950 
Critical Mach Numbers for Swept-Back Wings 


General Performance Reduction [Equations for 
Propellers 


The Flow Near the 
Centre of an Expansion 


Note on the Dependence of Flap Hinge Moment 
Derivatives on Hinge Position 4 


On the Chordwise Lift Distribution at the Centre of 
Swept Wings 


Erratum—(Volume I, Part vy) 
“A Note on the Bending Moment Induced in the 
Booms of a Spar at the — of cs of a 
Concentrated Load” ... 


PART III, NOVEMBER 1950 
Beams on Elastic Supports and on Cross Girders 
The Estimation of Range of Jet-Propelled Aircraft 
On Non-Steady Motion of Slender Bodies 


The Correction of Wind Tunnel Nozzles for Two- 
Dimensional, Supersonic Flow 


Some Developments of Methods for Solving 
the Flutter Equations 


Chang O’Chou 
J. Morris 


E. C. Pike 
John W. Miles 
Geoffrey Sewell 


S. Goldstein and 
G. N. Ward 


S. Neumark 


W.J. D. Annand 
N. H. Johannesen 
and R. E. Meyer 
W. J. Duncan 


D. Kiichemann 
and J. Weber 


H. F. Winny 


J. M. Klitchieff 
H. Pearson 
J. W. Miles 


R. E. Meyer and 
M. Holt 


J. Williams 


Printed by the Lewes Press (Wightman & Co. Ltd.), Friars Walk, Lewes, England, and published by 
Tue Royal AERONAUTICAL Soctety, 4 Hamilton Place, London, W.1, England. 


by 
eee eee 
eee eee 
eee 


